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INTRODUCTION 


This Annual Progress Report contains investigations on detonation 
driven magnetogasdynamic and solid state power generators and related 
research problems. The work was conducted under ONR Contract N00014- 
81-WR-10107 in the period from February 1 to September 30, 1981. 

The electric currents and voltages produced by electromagnetic 
induction in explosion driven magnetogasdynamic generators with homo¬ 
geneous external magnetic fields B q are calculated for two shock flow 
models, (1) the (implosion produced) jet flow v(t) » V Q H(t) with con¬ 
stant shock speed v q , and (2) the plane detonation flow with shock 
-1/3 

speed v(t) « t (similarity solution). External load circuits with 


resistance R and Inductance L are considered which are connected 
o o 

i) to the upstream and ii) to the downstream ends of the generator 
electrodes so that the magnetic self-fields B of the electric currents 


are i) antiparallel and ii) parallel to the external 


-ic field 


1 Q . It is shown that optimum power is generated for generators with 
jet shock flow and load circuits of type (ii), i.e., positive super¬ 
position of induced and external magnetic fields. In this case, elec¬ 
tric energies AE - 10 J Joule can be generated with external fields / 
B q » 1 Tesla per detonation. The importance of the results in t iK* N 
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Chapters II and III for the design of detonation driven generators is 
obvious. 


One of the main technical problems of explosion driven magneto- 
g a adynamic generators is the requirement of large external magnetic 
fields B q > 1 Tesla, which are difficult to procure in connection with 
military applications. For this reason, a solid state power generator 
is analyzed in which the driving electromotoric force is due to electric 
polarization in the stress field of an explosion produced shock wave. 

The solid is assumed to be a one-dimensional slab between plane elec¬ 
trodes , which are connected through an external load circuit. The tem¬ 
poral current build-up in the external circuit, due to electric polari¬ 
zation and electric conduction behind the shock front of the stress 
wave, space charge and electric polarization relaxations, is determined 
by an inhomogeneous Volterra integral equation. For shocked solids with 
cross section of the order of one square meter, energies of the order 
AE * 10^ Joule per pulse can be achieved. This scheme for an electric 
generator, which does not require external magnetic fields, appears to 
be promising for experimental Investigation (Chapter IV). 

An other method of avoiding large external magnetic fields B q > 

1 Tesla is to build an explosion driven MGD generator, in which a weak 

-2 

external magnetic field B q ~ 10 Tesla is compressed between the plasma 

shock front (high magnetic Reynolds numbers R ■ u ova) and a copper 

in o 

liner. The exact relativistic theory of the flux compression MGD gen¬ 
erator leads to two coupled integrodifferential equations which appear 
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to be unsolvable by analytical methods. For this reason, the electromagne¬ 
tic diffusion and induction in accelerated conductors in external li fields 

o 

was analyzed in the parabolic diffusion approximation, which is valid for 

"18 7 "1 

small relaxation times t ■ t Q /a (t ~ 10 sec for copper, a m 6x10 0 /m). 

The conventional boundary conditions (which violate the continuity of the 
tangential electric field ) were replaced by new boundary conditions, which 
consider the electromagnetic transients outside of the conductor. Thus, a 
realistic evaluation of magnetic flux compression and magnetic flux losses 
through the plasma piston and the copper liner is accomplished (Chapters V 
and VI). 

In Chapter VII, a quantum-kinetic theory of the anomalous plasma 
heating by high-frequency electromagnetic fields is presented based on the 
Maxwell-Vlasov and Schroedinger equations. This Investigation has been 
carried through in collaboration with visiting Prof. S.H. Kim. 

The work on flux compression electric power generators will be com¬ 
municated in the 1982 Annual Report. There, also new types of explosion 
driven generators will be discussed, which do not require an external 
magnetic field. 
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DETONATION DRIVEN INDUCTION GENERATOR 
WITH ANTIPARALLEL EXTERNAL AND INDUCED MAGNETIC FIELDS 


H.E. Wilhelm 

Michelson Laboratory, Naval Weapons Center, 
China Lake, CA 93555 


Based on Maxwell's equations, the electric current and voltage 
pulses induced in explosion driven induction generators * u Q avL >> 1) 
with plane electrodes and homogeneous external magnetic field B q are cal¬ 
culated analytically for two plasma shock flow models, (i) the jet flow 

with shock speed v(t) ■ v Q H(t) and (ii) the plane detonation flow with 

1/3 -1/3 

shock speed v(t) * (2/3)(E 0 /p Q ) t . The external load circuit with 

resistance R and inductance L is connected to the flow entrance ends of 
o o 

the electrodes so that the magnetic self-field B(t) of the generated cur¬ 
rent I(t) is antiparallel to the transverse external field B Q . It is shown 
that the jet flow produces large current and voltage pulses with extended 
plateaus, whereas the detonation flow generates considerably smaller and 
shorter current and voltage pulses. In both cases, the magnetic self-field 
B(t) of the generated current I(t) is of the order-of-magnltude of ^ o , and 
its negative superposition on the external magnetic field is not advanta¬ 
geous for power generation. 
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INTRODUCTION 


In contrast to the conductive magnetohydrodynamic generator with steady-state 
magnetic and flow fields —^, induction generators are based on the electromagnetic 

interaction of a conducting flow field with a magnetic field, where at least one of 
2 ) 

the fields is transient — so that voltages are induced due to temporal flux changes 
d$/dt J 0. In some induction generators, a detonation plasma with a decreasing shock 
velocity v(t) -*■ 0 for t -*■ • is directly injected into the inter-electrode space — , 
whereas in more recent experiments^*^discontinuous plasma jets are shot across an ex¬ 
ternal, constant magnetic field. The plasma jets v(t) ” v Q H(t) of nearly constant 
shock speed v q are produced, e.g., by Implosion of a dense noble gas between flyer 
plates and subsequent expansion through a diaphragm into the generator channel —\ or 
by means of a tube like arrangement of the solid explosive resulting in a cylindrical 
implosion with large axial shock velocity of the filling gas —K 

l 5 

Detonation driven induction generators have shock velocities v ■ 10 - 10 m/sec, 

4 -1 

external magnetic fields B q ■ 1 - 10 Tesla, conductivities a > 10 £1 /m, and flow 

1 2 10 12 _3 a n ^ 3 

energies pv '10 - 10 Joule m at pressures p > 10 Dar (p > 10 xg m ). The 

corresponding magnetic Reynolds number R^ is large, whereas the magnetic interaction 
number R^ is small. 


R_ * u ovL » 1 
ffl o 


R - B 2 /v pv 2 << 1 
1 o o 


where L * [Vx?/B|. For R^ » 1» the electric current is restricted to a thin layer 

6 - (t/pa) behind the shock front (t ■ transit time of flow in B ). In view of 

o o 

1 2 2 6 8 
the large dynamic pressure j pv compared to the magnetic pressure B Q /2u o * 10 - 10 

-3 

Joule m , the efficiency of energy transformation is small, while the performance is 

A 

significant since currents I - (B o /u o )b > 10 Amps (b - electrode width) and energies 

__» 4 4 5) 

Ult > 10 Joule have been produced experimentally per detonation —‘'. 
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Herein, the Initial-value problem for a detonation driven generator 

is solved in closed form by means of Maxwell's equations, in which the flow 

entrance ends of the coplanar electrodes are connected by an external L -R 

o o 

load circuit (Fig. 1) so that the magnetic field B(t) of the generated cur¬ 
rents I(t) is antiparallel to the external magnetic field B q . Two types of 
plasma shock flows are considered corresponding to i) the idealized jet 
flow with constant shock speed v q and 11) the plane detonation flow due to 
an explosive energy release E q per unit area in a gas of initial density p q : 

v(t) - vH(t) , x(t) - v tH(t) , 

o o 

0 < t < t , t - x co /v Q ; 

and 

v(t) - (2/3)(E o /p o ) 1/3 t _1/3 , x(t) - (E o /p o ) 1/3 t 2/3 

0 < t < t , t - (E /p )' 1/2 x 3/2 

o o 00 

where x(t) is the shock front position, and t is the transit time of the 
shock front x(t) in the duct of effective length x v (Fig. 1). The self¬ 
similar solution (3) is asymptotically correct for larger times t —'', but 
can be used for t » 0 since the singularity of v(t) at t * 0 is without 
physical consequences for generator applications. 

The external circuit is assumed to be located at x ■ -x (Fig. 1), 

o 

where x is an effective distance which can be calculated once the geometries 
o 

of the inductance L and resistance R are specified, 
o o 


( 2 ) 


(3) 
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ieid B q , and R q - L load circuit AD. Shock front 






INITIAL-VALUE PROBLEM 


Fig. 1 exhibits a one-dimensional model for a detonation driven generator, in 

which the shock wave flow v(t) * dSt(t)/dt is initiated at time t = 0 in the plane 

x ■ 0 by explosive charges in the space -x < x < 0 so that the shock front is at 

x(t) > 0 for t > 0. The Cu-electrodes in the planes z = ±c/2 have the extensions 

-x q < x < x^ < ® and -b/2 < y < +b/2, and are connected by an external circuit AD 

with load R and inductance L . As a result of the vxli interaction of the conducting 
o o o 

shock front (o > 10^0 Vm) with the external magnetic field ? * {0 ,B q , 0}, a current 

I(t) is induced behind the shock front BC in a cross section AxAy - b6 where 
A 1/2 

6 - (t/p Q o) is the (small) width of the current layer. The current I(t) flows 
in the closed loop A+ B+ C* D + A, where the leg BC at R(t) moves with the shock 
speed v(t) * dx(t)/dt. 

Application of the surface curl equation nx[5] • p Q t(t)/b to the interfaces AB, 
BC, CD, DA demonstrates - in the one-dimensional approximation neglecting end effects 
at y * ±b/2 and x * -x q - that the current I(t) produces a magnetic self-field in the 
extending loop volume, 

f(t) - {0,-u Q I(t)/b,0} , -x Q < x < x(t) , |y| < y, M < f t W 

which is antiparallel to the external magnetic field The total magnetic field 

within the loop volume is 

!(t) - {0,B o -v o I(t)/b,0} , -x Q < x < x(t) , |y| < | , |z| < f . (5) 

whereas $(t) » t outside of this region. Integration of the Maxwell equation Vxi ■ 
o 

-al/3t over the cross secti' ABCD of le extending current loop yields by Stokes' 
theorem 

$l*-dr - - 0f| + <P <$**)•& 
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(6) 




where 


v << 1/p e 

o o 


(7) 


and v(t) is the velocity field of the sections of the closed loop, at the point 
corresponding to the path element dr[v(t) ■ e x dx(t)/dt for BC, but v(t) * t) for CD, 
DA, AB]. The surface element dS * -e^cdx forms with dr [pointing everywhere in the 
direction ^(t)] a right hand system. By Eqs. (6) and (7) 


(R o + Rp)I(t) + L 


iliEl « _ c J (t> + J c/Z [B + B(t)]dz 

-x dt -c/2 dt ° 


'o dt 


( 8 ) 


since / E *dz * IRp (R„ * plasma resistance), / Edz * R q I + L Q dI/dt, and /Edz * 0 

® D s —l 

for the sections AB and CD of the highly conducting electrodes (o Cu - 10 £5 /m). 

Accordingly, 


where 


(R o + yi(t) + £ { [L o + i.(t)JI(t)} - cB o ^i^- 


L(t) * p c[x + x(t)]/b >0 , 0 < x(t) < X 

o o 


(9) 


( 10 ) 


is the time-dependent inductance of the electrode sections AB & DC of increasing 
length Ax - x Q + x(t). Multiplication of Eq. (9) by I(t) leads to the conservation 
law 


<R 


+ v i2 + [L - 


dt 2 


+ L(t)]i 2 } + ^i 2 -i2M 




(ii) 


according to which the ohmic power dissipation, the temporal change of the magnetic 
energy, and the work on the current sheath per unit time equal the power input 01. 
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The plasma resistance Rp * Rp(t) is a slowly varying function of t and can be 
calculated only via a detailed physical model of the current sheath 6 at x % x(t). 
However, since R q is sufficiently large in actual experiments, R q + Rp can be treated 
as a quasi-constant, 

R q + Rp(t) * constant, R q > Rp ; R Q + Rp(t) - R q , R q » Rp . . (12) 

For a generalized formulation, dimensionless dependent and independent vari¬ 
ables are introduced by 


J(t) - I(t)/I , t = t/t , £(t) ” x(t)/c 


*0 * - <v 2 /b)/«o + V . » 0 - V c > + V ( V 2/b) - (U) 

By Eqs. (9) - (10) and (13) - (14), the electric current I(t) in the explosion driven 
generator is determined by the inhomogeneous initial-value problem with variable co¬ 
efficients: 


{[A o + £(t)]J(t)} + J(t) - , 0 < t < t 


J (t - 0) - 0 


, (15) 


Eqs. (15) - (16) are solved for the jet and detonation flows in Eqs. (2) and 

(3), respectively, which define the dimensionless coordinate CCO of the shock front. 

12 2 

In view of the large flow energy, y pv » B q /2p o , it is not necessary to correct 
the above generator model for velocity decreases due to Lorentz and viscous forces. 
Furthermore, viscous boundary layers are disregarded since the plasma moves in form 
of a rectangular velocity slug through the duct in large Reynolds number flow, 

R ■ pvL/y >> 1. 
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JET-FLOW SOLUTION 


For the jet flow, the dimensionless shock front coordinate £(x), C Q » and the 
transit time t are by Eqs. (2) and (13) 


v p c /b 

£(t) - E q t , ) 


x R+L 
-2L(-Q—L.) 
v ' 2.J 

O U Q C IQ 


(17) 


Since ECO > 0 for t > 0 , the initial-value problem (15) - (16) is without singu¬ 
larities, and has the closed-form solution: 


e e -(i+E rt )/E„ 

j (t) * x?f~[i " (i + 5T t) ] » o<t<; 


(18) 


The transit time in Eq. (17) satisfies the inequality, 

L 


(E /A )t * ^—1 >:> 1 > for x » X » L «C (x /c)p c 2 /b , (19) 


p c 2 /b 


which indicates that the independent variable (£ q /A q )i may assume large values. 
Eq. (18) has the limits: 


J(0 - (E /A )x , (E /A )t << 1 

O O O 0 


J(t) - E„/(l + EJ , (E M )r » 1 


( 20 ) 

( 21 ) 


Accordingly, the maximum current is reached only asymptotically for large times, 
since J(t) is a monotonically increasing function of 0 < t < t, 

B 


max o 


€«/<! +O » I - ^ b^/d + E ), (E/A)t»1 . (22) 


max p 
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The generated voltage, U * R q I + L o dI/dt, is by Eq. (18) in dimensionless representa¬ 
tion 


5 -(1+E )/E L /R £ 5 0 -(1+25 )/5„ 

u(t) a + ^T) 

o o o o 


where U • I R . Eq. (23) has the asymptotic values 
o o o 

L /R 1 + 25 5 L /R 5 5 

« i , (24) 

o o o o o o 


"< T >-rrr * r T>>1 

o o 


For 1 + (5 /A )t - (L /R t )(1 + 25 )/A , U(r) assumes a minimum, 
o o o o o o o 


5 L /R 1 + 25 -d+5 )/£ 

o r o o , o, 1 00 


min 1+5 


-2— _Sd_ r- u — (-Si) ] 

+ 5 1 + 25. 1 t V A. 


Accordingly, the maximum voltage value is reached only asymptotically for large times, 

(5 fl /A o )T >9, 


U -5/(1 +5), U «IR 5„/(l + 5 ) 
max o o max o o o o 


Example . For b * 10 ^m, c ■ 1 m, x ■ 10 ^m, « 1 m, a • 10 3 0 ^/m, Rp - 

o -1 c/b6 - 10 -2 0, R + R_ - 2xl0 _1 fJ, L - 4»«10~ 6 Henry, v - 10*m/sec, B - 1 Tesla, 
o T o o o 

-5 «. -4 

the characteristics of the generator are: t Q » 2irxl0 sec, t - 10 sec, 5 Q * 2n/10, 

A - 11/10, I - 10 6 /4v Amps, I 1 3x10* Amps, U_- 6xl0 3 Volt, (UI) „ i 2x10® Watt, 

o o max max 

_„ 4 

AE - Ult - 10 Joule. 
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Fig. 2 shows the current J(t) versus (£ /A )t for £ ■ 2ir/10 and £ ■ 9/10 

o o o o 

based on Eq. (19). The distribution J(x) rises steeply within 0 < (£ o /A q )t < 1 and 
then flattens out into a shightly increasing "plateau" for larger times t. In ex- 

4 5) 

periments — 1 —, the 'plateau" is more constant or even decreases somewhat with in¬ 
creasing time t, which appears to be due to i) shock speeds v(t) decreasing slightly 
with time t, and ii) external circuit leads to the center of the electrodes. 







Current transient J(x) of jet shock flow versus (£ /A )t for 









DETONATION-FLOW SOLUTION 


The dimensionless shock front coordinate £(x), and the transit time x for 
the detonation flow follow from Eqs. (3) and (13) as 

» . 2/3 . , E o,l/3 -1, y o C /b .2/3 - .-3/2 , x ® x 3/2 

5(t)-5 0 t , 5 0 "(—> ‘ * T " 5 o ( T } * (28) 

Since £(x) > 0 for x > 0, the solution of the initial-value problem (15) - (16) is 
representable in terms of an indefinite Integral with regular integrand: 


J(t) - 1 - 


i + Cl, 


T 3C* 1 [ t 1 ^ 3 -a o 1 8Ectan(a 0 T 1 ^ 3 )I 


o J 0 


dr 


. 2 2/3. W; 1 tt 1/3 -°; 1 ar«a»C V l/3 )] 

(1 + o q t )e 


where 


a o S W 


1/2 


, 0 < t < r , (29) 


, (30) 


1/3 

Eq. (29) is an integral functional of a Q x , where 


V 1/3 m { T /l ~t + — 2~ ]}1/2 ?> 1 * for x . >> x o • L o 5 (* 0 /c)y 0 c 2 /b. 
pc / b 

° (31) 

For small and large times, Eq. (29) gives the limits 


• 1(0 


J(t) ■ 0 , a t >> 1 


1/3 «1 

9 

(32) 

1/3 »1 


(33) 


J(x) has a true m ax im um J at time x - x where dJ(x)/dx * 0 for x ■ f. By Eqs. (15). 

max 

(28), and (29), and ? are given by 
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(34) 


j - c /(c + 4 * 1/3 ) 

max o o 2 


2 _ , -1/3 

3 C o + T 


- , -l r 1/3 -1 . , 1/3,, 

, t 3? [t -a arctan(a r )] 

1 + A 1 / e ° ° ° dx 

_2_iQ_ 

,-l r -l/3 -1 , -1/3,, 

„ 2-2/3, 35 o [T -°o arctan( V )] 

(1 r a 0 T )e 


The generated voltage, U Q ■ R q I + L o dI/dt, is by Eq. (29) in dimensionless 
representation 

U(t) - J(t) - -^2 [(1 + 1 5 o t“ 1 / 3 )j(t) - -| 5 0 t“ 1 / 3 l//V 0 (l + a 2 x 2/3 ), 0 < t < t. 


The inductive singularity, U(x) -*■ » for r -*■ 0, is due to the divergence of the 
selfsimilar speed —^, d£(x)/dx - (25 q /3)t ^ 73 ® for x -*■ 0, and has, therefore, 

no physical meaning. 

Example . For b * 10 'Sn, c * 1 », x q » 10 ^m, x^ » 1 a, a ■ 10 3 ft Vm, 

Rp i a _1 c/b6 « 10~ 2 n, R o + Rp - 2xl0*’ 1 n, L o * 4irxl0~ 6 Henry, v(10 _6 sec) - 

4 -5 

2.650x10 m/sec, B ■ ITesla, the generator characteristics are: t * 2irxl0 sec, 
o o 

t - 1.262xl0~ 4 sec, C - 2ir/10, A * 11/10, I - 10 6 /4ir Amps, I * 1x10* Amps, 
o o o max 

U * 2x10 3 Volt, (01) * 2x10 7 Watt, AE - Ult ' 10 3 Joule, 

max max 

1/3 

Fig. 3 exhibits the current J(t) versus a Q x for 5 q * 2ir/10 with A q ■ l/4ir, 
10/4 tt, and 100/4x as parameter based on Eq. (29). In each case, J(t) rises rap¬ 
idly to a maximum value and then decays slowly with increasing time. The peak 

current J decreases with increasing A (inductance). It should be noted that 
max o 

»l/3 1/2 

the distributions J(x) are applicable up to the transit value o q x • (x^/A^c) 

[Eq. (28)] which depends on the respective system parameters. Experimental 

3) 

current oscillograms for the detonation flow— agree only qualitatively with 
Eq. (29) since the electrodes were connected in their centers by the external 


load circuit. 
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Current transient J(t) of detonation shock flow versus 
1/2 1/3 

(C 0 /A q ) t for £ - 2 tt/ 10, and A - 1/Air; 10/4*; 100/4* 










CONCLUSION 


An elementary, one-dimensional theory for an explosion driven in¬ 
duction generator has been developed, in which the upstream ends of the 

electrodes are connected by an external R -L load circuit. The formu- 

o o 

las derived for the current I(t) and voltage U(t) permit an optimum gen¬ 
erator design by appropriate selection of the system geometry, flow, and 
load parameters. 

The jet flow with (about) constant shock speed v(t) * v q produces 

an order-of-magnitude more energy than a comparable detonation flow with 
-1/3 

shock speed v(t) « t decreasing with time. Since shock speeds in¬ 
creasing with time are not realizable in the coplanar duct, the jet flow 
is ideal for power generation. 

The magnetic self-field B(t) of the generated current I(t) is of 
the same order-of-magnitude as the external magnetic field B q . B(t) 
considerably reduces the generator performance if the electrodes are 
connected at their upstream ends by the external load circuit (5(t) and 
1$ o antiparallel). 

An optimum power output is achieved by attaching the external load 

circuit to the downstream ends of the electrodes so that B(t) and B are 

o 

in the same direction. For a maximum energy output the parameters £ 0 > 

I Q , t should be chosen as large and A^as small as possible. 

For the transit period t of the shockwave, explosion driven gener¬ 
ators are generally contained by inertial confinement. For this reason, 
no attempt has been made to discuss generator performance for times 
beyond t. 
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DETONATION-DRIVEN INDUCTION GENERATOR 
WITH POSITIVE SUPERPOSITION OF EXTERNAL AND INDUCED MAGNETIC FIELDS 

H. E. Wilhelm 

Michelsoti Laboratory, Physics Division 
Naval Weapons Center, China Lake, California 93555 


ABSTRACT 

Based on Maxwell's equations, the electric current and voltage pulses 

induced in explosion-driven induction generators (R * p ovL >> 1) with plane 

mo 

electrodes and homogeneous external magnetic field B q are calculated in closed 

form for an external load circuit with resistance R and inductance L con— 

o o 

connected to the downstream ends of the electrodes so that the magnetic self¬ 
field B(t) of the generated current I(t) is in the same direction as the 
transverse external field B q . Two plasma shock flow models are considered: 

(i) the jet flow with shock speed v(t) - V Q H(t) and (ii) the plane detonation 

flow with shock speed v(t) »■ (2/3) (E /p )^^t It is shown that the jet 

o o 

flow produces large cumulative current and voltage pulses which increase with 
time, whereas the detonation flow generates considerably smaller and shorter 
current and voltage pulses. In both cases, the magnetic self-field B(t) of the 
generated current I(t) is of the order of magnitude of B q . In case of positive 
superposition of the Induced magnetic field on the external magnetic field B q 
maximum flux changes are obtained. For this reason, in induction generators 
the external load circuit should not be attached to the upstream ends or the 
centers, but to the downstream ends of the electrodes. 
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INTRODUCTION 


Induction generators are based on the electromagnetic interaction of an 
electrically conducting flow field with a magnetic field, where at least one of 
the fields is transient so that voltages are induced in an external circuit 
which are due to temporal magnetic flux changes d$/dt ^ 0. An explosion-driven 
Induction generator with a transient plasma shock flow in a constant external 
magnetic field B q has been analyzed^ when the external R q - L q load circuit is 
attached to the flow entrance ends of the electrodes. In this case, the self- 
magnetic field B(t) of the generated current I(t) is antiparallel to the trans¬ 
verse external magnetic field B q so that the magnetic flux changes d$/dt during 
the transit time t of the flow are reduced.—^ Herein, a novel explosion-driven 
induction generator with optimum magnetic flux changes is proposed, in which 

the external R - L load circuit is connected to the flow exit ends of the 
o o 

coplanar electrodes so that the Induced magnetic field B(t) of the current I(t) 

in the shock front current layer and the electrodes can superimpose themselves 

positively on the external magnetic field B q . This effect is quantitatively 

significant since B(t) and B q are parallel and of the same order of magnitude. 

Thus, considerably larger current, voltage, and power pulses are induced than 
-*■ -+■ 

in the antiparallel B q - B case. 

In induction generator experiments, the external R q - L q load circuit is 

2 3 4) 

usually attached to the center of the electrodes.— 1 —In earlier experiments, 

a detonation-produced plasma with a decreasing shock velocity v(t) -*■ 0 for t -*■ » 

2 ) 

is directly Injected into the inter-electrode space,— whereas, in more recent 
3.4) 

experiments,— 1 — discontinuous plasma jets are driven across an external, 
constant magnetic field B q . Plasma jets v(t) * v Q H(t) of nearly constant shock 
speed v Q are produced, e.g., by implosion of a dense noble gas between flyer 
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3) 

plates and subsequent expansion through a diaphragm into the generator channel,— 
or by means of a tube-like arrangement of the solid explosive resulting in a 
cylindrical implosion with large axial shock velocity of the filling gas.—^ 

Detonation-driven induction generators have been operated with shock veloci- 
4 5 

ties v * 10 - 10 m/sec, external magnetic fields B * 1 - 10 Tesla, conductivi- 

o 

4 “1 1 2 10 12 —3 

ties a > 10 0 /m, and flow energies ^ pv ' 10 - 10 Joule m at pressures 

4 2-3 

p > 10 bar (p > 10 kg m ). The corresponding magnetic Reynolds number R 

m 

is large, whereas the magnetic interaction number R^ is small: 

R = u ovL >> 1 , R. - B^/p pv^ << 1 , (1) 

mo loo 

where L ' Ivxb/bI. For R » 1, the electric current is restricted to a thin 
1 m 

A i/2 A 

layer 6 ~ (t/p Q a) behind the shock front (t ■ transit time of flow in B q ). In 

1 2 

view of the large dynamic pressure ^ pv compared to the magnetic pressure 
2 6 8 —3 

B*/2p q ■ 10 - 10 Joule m~ , the efficiency of energy transformation is small, 

while the absolute performance is considerable since currents I * (B /p )b > 

o o 

4 —* 4 

10 amps (b ■ duct width) and energies Ult > 10 Joule have been obtained experi- 

3,4) 

mentally per detonation.— 2 — 

In the following, the initial-value problem for a detonation-driven generator 
is solved in closed form by means of Maxwell's equations, in which the downstream 
ends (x * x q ) of the coplanar electrodes are connected by an external L q - R q 
load circuit. In this case, solutions for the magnetic field B(t) of the 
generated currents I(t) exist which are parallel to the external magnetic field 
B q (Fig. 1). Two types of plasma shock flows are considered corresponding to 
1) the Idealized jet flow with constant shock speed v q and ii) the plane detona¬ 
tion flow due to an explosive energy release E q per unit area in a gas of 
initial density p Q : 

v(t) - v Q H(t) , x(t) - v Q tH(t) , (2) 

0 < t < t , t » x Q / v 0 * 
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load circuit CD. Shock front AB at 





and 


v(t) = (2/3)(E o /p o ) 1/3 t" 1/3 , x(t) = (E o /p Q ) 1/3 t 2/3 , (3) 

0 < t < t , t = (E /p ) _1/2 x 3/2 

o o o 

where x(t) is the shock front position, and t is the transit time of the shock 
front x(t) in the duct of effective length x q . The self-similar solution (3) 
is asymptotically correct for larger times t,—^ but can be used for t > 0 since 
the singularity of v(t) at t * 0 is without physical consequences for generator 
applications. 

The mechanical stability of the generator during the shock transit time t 
is provided for by inertial confinement of the duct through high density plaster. 
No attempt is made at discussing the electrical output for times t > t of the 
explosive disintegration of the generator. The external circuit location x q is 
an effective one, which can be calculated once the geometries of the inductance 
L and resistance R are specified. 
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INITIAL-VALUE PROBLEM 


Figure 1 exhibits a one-dimensional model for a detonation-driven generator 

in which the shock wave flow v(t) = dx(t)/dt is initiated at time t = 0 in the 

plane x = 0 by explosive charges in the space x < 0 so that the shock front is 

at x » i(t) > 0 for t > 0. The copper electrodes in the planes z = ±c/2 have 

the extensions x > x q < 00 and -b/2 < y < +b/2 and are connected by an external 

circuit CD with load R and inductance L at the flow exit end x = x . The 
o o o 

entire system is embedded in a homogeneous magnetic field B q * {0,B o >0}. As a 
result of the v x B q interaction of the conducting shock front (a > 10^ Q Vm) 

with the external magnetic field B q , a current I(t) is induced behind the shock 

* 1/2 

front AB in a cross section AxAy = b6, where a ~ (t/u Q o) is the (small) width 

of the current layer. The current I(t) flows in the closed loop A -+• B -*■ C -► 


D -*■ A, where the leg AB at x * x(t) moves with the shock speed v(t) * dx(t)/dt. 


Application of the outface cutl equation n x fS] - » 0 t(t)/b to the inter- 
faces AB, BC, CD, DA demonstrates — in the one-dimen3ional approximation 
neglecting end effects at y = ±b/2 and x * x q — that the current I(t) producer 
a magnetic self-field in the extending loop volume. 


B(t) * {0,u o I(t)/b,0} , x(t) < x < x Q , |y| < , | z| < , (4) 

which is parallel to the external magnetic field if I(t) > 0. The total 
magnetic field within the loop volume is the positive superposition 


B(t) - {0,B o +u o I(t)/b,0} , x(t) < x < x Q , |y{ < j , [z| < j , (5) 

whereas ^(t) ■ B q outside of this region. Integration of the Maxwell equation 
V x | - -3l/3t over the cross section ABCD of the extending current loop yields 
by Stokes' theorem 
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•dr * - $ • dS + <p (v x b) *dr , 


( 6 ) 


where 


E * E + v x B 


v << 1/y e 
o o 


(7) 


and v(t) is the velocity field of the sections of the closed loop, at the point 
corresponding to the path element dr [v(t) = e x dx(t)/dt for AB, but v(t) » 0 for 
CD, DA, CB]. The surface element d^ = +e^cdx forms with dr [pointing everywhere 
in the direction ?(t)], a right-hand system. By Eqs. (6) and (7) 

«„ + V I(t > + L „ ^ - -= l ° :+ + »<«w* • < 8 > 

x(t) -c/2 

B * D 

since / E dz * IRp (Rp = plasma resistance), / Edz ■ R q I + L Q dI/dt, and /Edz 4 0 
for the sections BC and DA of the highly conducting electrodes (a ~ 10 ft /m). 
Accordingly, 

(R q + Rp)I(t) + {[L o + L(t)]I(t)} - cB q , 


(9) 


where 


L(t) * p c[x - x(t)]/b > 0 


0 < x(t) < x 


( 10 ) 


is the time-dependent inductance of the electrode sections AD and BC of decreasing 
length Ax ■ x q - x(t). Multiplication of Eq. (9) by I(t) leads to the energy 
conservation law 


(E o + R ,,)! 2 + i {i (L 0 + £(t)]I 2 } 


o _2 dx(t) _ dx(t) _ 

2b“ 1 ~d r~ + 1 ~dt B o c 


( 11 ) 


The plasma resistance Rp *• Rp(t) is a slowly varying function of t and can 
be calculated only via a detailed physical model of the current sheath 5 at 
x * x(t). However, since R q is sufficiently large in actual experiments, 

R q + Rp can be treated as a quasi-constant. 
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R q + Rp(t) a constant , R q > Rp ; R q + Rp(t) = R q , R q >> Rp . I 

For a generalized formulation, dimensionless dependent and independent 
variables are introduced by 


J(t) = I(t)/I Q , x = t/t o , C(t) = x(t)/c 


• < B o / “o> b • 'o ■ < u o c2/b)/<R o + V • A o * <V c) + ‘■o /( “o c2/b) • (14> 

By Eqs. (9), (10), (13), and (14), the electric current I(t) in the explosion- 
driven generator with positive external and induced magnetic field superposition 
is determined by the inhomogeneous initial-value problem with variable coeffi- 


{[A - £(t)]J(t)} + J(t) 


, 0 < t < t 


J(r - 0) - 0 . (16) 

Equations (15) and (16) are solved for the jet and detonation flows given by 
Eqs. (2) and (3), respectively, which define the dimensionless coordinate |(t) 
of the shock front. 

12 2 

In view of the large flow energy, -j pv >> B q /2vi o , velocity decreases due 
to Lorentz and viscous forces are neglected in the above generator model. Further¬ 
more, viscous boundary layers are disregarded since the plasma moves in the form 
of a rectangular velocity slug through the duct in large Reynolds number flow, 

R ■ pvL/u » 1. 








3T 


T 


JET-FLOW SOLUTION 


For the jet flow, the velocity and coordinate of the shock front are given 

in Eq. (2). The dimensionless shock front coordinate £(t), 5 , and the transit 

o 

time t are, by Eqs.(13) and (2), 


x- R- + Rp 


v u c / b 

* «o-T<rnr> * 

O r O U c /b 


) . 


(17) 


I 

I 

I 

I 


The initial-value problem (15) and (16) has for C(t) - £ t a closed-form solution. 

o 1 

which is well behaved: 

J(t) - y~ ' ° f [1 - (1 - jp t) ( 1 ”V / 5 o] , 0 < t < t , (18) 


where 


lim J(t) - - in(l — t) < 

- 1 


T < 1 


(19) 


The critical time t defined by £ t/A - 1 is not attainable within the transit 

00 

period t since 


Ti r M(1 + 7-ex/ b 


■)t , t < t for L > 0 


( 20 ) 


For short times t -*■ 0 and the maximum time t, Eq. (18) gives, respectively, 
the limits 


J(x) - a o /A Q )T , (C 0 /A 0 )t « 1 , 

o o 


( 21 ) 

( 22 ) 


Equation (22) represents the maximum generator current for arbitrary £ q > 0 and 
leads to the special formulas 
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J max S V (1 ‘ V for C o < 1 » L o << W o CX o /b * 

J max ’ “ * n(1 - r S> < - f or E - 1 . 


(23) 

(24) 


The generated voltage, U = R^I + L^dl/dt, is by Eq. (18) in dimensionless repre¬ 
sentation (U = i R ) 
o o o 

u<t> ■ r^r 11 - (1 -r ’> (1 ' E ° )/5 °i ♦ «• - , 

o o o o o 


o < X < T , (25) 


where 


L /R 

lim U( t) - - *n(l - t) + °— A* 1 (1 - A~ t) -1 < ®, 0 < t < t 

5 -*■ 1 o o 0 0 


(26) 


and 


U( x) * [1 +-2—2- ( 


L„/R„ 1 - 2C 5 L /R E E 

— )] 7^ t +-2—. 1 r-T«l , (27) 


t A ’ A 
o o o 


U 


C ° ri n C ° , L ° /R ° n ^ (1-2EO/E,., 

^ g 11 - (1 - — t) o oj + —-— (1 - — T ) O o , 

^o o o A o o 

T - T . ( 28 ) 

For special cases, the maximum voltage is given by simpler formulas, 

U m.x * V (1 -V fOC V 1 • L o ~ “o c V b ■ < 29 > 

“ mx ' - to < 1 - C" + t A o 1<l - A ; lf )' 1 * * £ “ o " 1 ■ (30 > 


Example. For b ■ 10 ^ m, c ■ 1 m, x q ■ 1 m, o * 10"* £2 Vm, Rp a o ^"c/b6 - 

10~ 2 £2, R q + Rp - 2 x 10 _1 £2, L q - 4ir x 10 _b Henry, v q - 10 4 m/sec, B q - 

1 Tesla, the characteristics of the generator are: t Q ■ 2ir x 10 ^ sec, t ■ 

10' 4 sec, - 2ir/10, A q - 2, I q - 10 6 /4ff amps, I nax : 1 * 10 5 amps, U mar * 

2 x 10 4 volts, (UI) * 2 x 10^ watts, AE ■ Ult - 10^ joules. 

max 
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2ir/10 and £ 


Figure 2 shows the current J(t) versus (£ o /A q )t for £ q - 

9/10 in the case of L << y cx /b so that t g t based on Eq. (18). While for 
o o o 

£ q = 2tt/ 10 the current rises steadily to J max * 2w/(10 - 2r) * 1.690, J(t) ~ 1 

nearly throughout the interval 0 < t < t for £ q * 9/10, and the sharp rise to 

J “9 occurs not until the very end t s t of the flow transit period. It 
max 

should be noted that the generator produces in each case average currents 

I > I in excess of the normalization current I ■ (B /u )b. Experimental 
o o o o 

current oscillograms are not available for comparison, since the advantages of 
a generator with a load circuit attached to the downstream ends of the elec¬ 
trodes have apparently not yet been recognized. 
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J 

max 


ViV*,/' 2 - £ o> fot £ < 


i ch 2 ' 3 * 1 ' 3 

4 O 


L « p cx /b 
o o o 


(37) 


Eq. (32) has positive, monotonically increasing solutions J(t) > 0 for z e(Q,t). 
The magnetic field 1 (t) of this current is in the same direction as the external 
magnetic field B q (positive superposition). 

The generated voltage, 0 ■ R Q I + L Q dI/dt, is by Eq. (32) in dimensionless 
representation (U q » I q R o )» 


L /R (2/3)£ - [t 

U( t) - J(t) + -S- 2 {- 2 - 


1/3 


. „ 2 2/3. 1/3 

A (1 - o t )t 


- (2/3)5 0 ]J(t) 

f 


0 < T «S T 


(38) 


The inductive singularity, U(r) -*■ « for t -► 0, is due to the divergence of the 
self-similar speed,—^ dCCO/dr - (2 £ q /3)t -*• • for t ■* 0, and has, therefore, 

no physical meaning. 

Example. For b » 10 ^ m, c ” 1 m, x * 1®, o » 10 3 ft 3 /m, R * o ^c/b6 * 

o p 

10 -2 fl, R q + R p - 2 x 10 _1 0, L q - 4ir x 10~ 6 Henry, v(t * 10~ 6 sec) - 2.650 x 

4 -5 

10 m/sec, B q ■ 1 Tesla, the generator characteristics are: t * 2ir x io sec, 

t * 1.262 x io 4 sec, £ * 2ir/10, A ■ 2, I * 10^/4r amps, I » 3xl0 4 amps, 

o o o max 

0 * 6xl0 3 volts, (UI) * 2xl0 8 watts, AE * Ult ' 10 4 joules, 

max max J 

1/2 1/3 

Figure 3 shows the current J(t)~ versus (5 /A ) t for £ * 2ir/10, 

o o o 

£ q ■ 9/10, and A q * 2 based on Eq. (32). In each case, J(t) increases monoto- 

» 1/2 1/3 

nically from 0 to its final value J at transit time t, where (£ /A ) t * 

max o o 

1 M 

(x /A c) • 1^2/2 for x ■ 1 m, c ■ 1 m, and A ■ 2. Cfa the average, the pro- 
o o o o 

duced currents 1 < I are smaller than the normalization current I » (B /u )b. 

o o o o 

Comparison with Fig. 2 indicates that J(t) of the detonation flow is by an order 
of magnitude smaller than J(t) of the jet flow. 
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CONCLUSIONS 


An electromagnetic theory for an explosion-driven induction generator has 

been developed in which the downstream ends of the electrodes are connected by 

an external R - L load circuit. The formulas derived for the current I(t) 
o o 

and voltage U(t) permit an optimum generator design by appropriate selection 
of the system geometry, flow, and load parameters. 

The jet flow with (about) constant shock speed v(t) » v q produces consid- 

_l/i 

erably more energy than a comparable detonation flow with shock speed v(t) <* t 

decreasing with time. Since shock speeds increasing with time are not realizable 

in the coplanar duct, the jet flow is ideal for power generation. 

The magnetic self-field B(t) of the generated current I(t) is of the same 

order of magnitude as the external magnetic field B . The induced magnetic 

o 

field 5(t) considerably increases the performance of generators with shock 
flow, since it superimposes itself positively on the external magnetic field B 
if the electrodes are connected at their downstream ends by the external load 
circuit (B(t) and B q in the same direction for optimum magnetic flux changes). 

An optimum power output is achieved by (i) attaching the external load 
circuit to the downstream ends of the electrodes and (ii) driving the generator 
with a jet shock flow. For a maximum energy output, the parameters £ 0 » I o , t 
should be chosen as large as possible, and A q should be chosen as small as 
possible. 
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POWER GENERATION BY SHOCKWAVE INTERACTION WITH SOLIDS 
ELECTRICALLY POLARIZABLE BY STRESS 

H.E. Wilhelm 

Michelson Laboratory, Naval Weapons Center, China Lake, CA 93555 

Abstract 

An elementary theory of the interaction of shock waves with 
solids which are electrically polarizable by stress is given and 
discussed with regard to the generation of electric power pulses. 
The solid is assumed to be a one-dimensional slab between plane 
electrodes, which are connected through an external load circuit. 

It is shown that the temporal current build-up in the external cir¬ 
cuit, due to electric polarization and electric conduction behind 
the shock front of the stress wave, space charge and electric po¬ 
larization relaxations, is determined by an inhomogeneous Volterra 
integral equation. The latter is solved analytically, and current 
pulse forms for incident step-shaped stress waves are calculated. 
For shocked solids with cross section of the order of one square 
meter, energies of the order AE - 10^ Joule per pulse should be 
achievable. 
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INTRODUCTION 

4 

Although electric energy pulses of 10 Joules and larger have been 

1-3) 

produced in explosion driven induction generators - , the requirement 

of large external magnetic fields of at least 1 Tesla renders these mag¬ 
neto gas dynamic generators unsuitable for many applications, in particular 
since the prospects for self-excitation of such magnetic fields during 
pulse times of the order 10 ^ 3ec are not promising . For these reasons, 
other methods for conversion of chemical energy of explosives into elec¬ 
tric energy have to be found. E.g., such a "non-magnetic" power supply 
could be used to produce the magnetic field required for magnetogasdyna- 
mic generators. 

Not only piezoelectric crystals such as tourmaline or ferroelectric 
crystals such as lead circonate titanate, but also a large number of 

other solids become strongly electrically polarized when subject to pres- 
4 

sures p > 10 bar. Such materials are, e.g., electrets, semiconductors, 

4) 

alkali halides, plexiglas, polysterene — . These solids simultaneously 
develop good to metallic electrical conductivities when subject to pres- 

4 6 

sure p ■ 10 -10 bar, due to nonthermal ionization by overlapping of at¬ 
omic wave functions at extreme pressures —^. Similar to the chemical 
polarization field in a battery, the electric polarization caused by the 
stress shock wave in the solid drives an electric current pulse through 
an external circuit attached via electrodes (transverse to the polariza¬ 
tion vector) to the solid. 

The interaction of a shock wave with solids, which are electrically 
polarizable by stress, appears to have potential for the generation of 
electric power pulses. Table I gives the stress polarization coefficient —^ 

f, the dielectric permittivity —^ e ■ e*e Q , the density p, and (dimensional 
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kJ±Jr. 




TABLE I: Order-of-magnitude estimates of stress polarization P, electric 

2 

field E, current I, and energy AE for a stress polarizable slab (A = 1 m , 

-2 3 

a =» 10 m) subject to a shock wave with speed v * 10 m/sec. 



Quartz 

BaTi0 3 

Rochelle Salt 

f[A sec/N] 

-12 
2.3x10 i 

l.OxlO -10 

4.0xl0~ 10 

p[kg m 3 ] 

2.65xl0 3 

6.02xl0 3 

1.77xl0 3 

a [N/m 2 ] 

| 1.3xl0 9 

3.1xl0 9 

0.9xl0 9 

[A sec/Vm] 

4x10 -11 

2xl0’ 8 

9xl0 _11 

2 

P[A sec/m ] 

3.0xl0~ 3 

3.1xl0 _1 

3-lxlO' 1 

E[V/m] 

7.6xl0 7 

1.5xl0 7 

3.9xl0 9 

I [A] 

3xl0 2 

3xl0 4 

3xl0 4 

AE[J1 

lxlO 5 

5xl0 6 

lxlO 8 
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„ x 2 

analysis) order-of-magnitude estimates of the stress amplitude a ~ ypv , 

electric polarization P ■ af, electric field E - P/e, short circuit current 

2 

I ' APv/a, and electric energy AE - (P /2e)Aa for one piezoelectric and two 

2 -2 

ferroelectric slabs of cross section A - 1 m and thickness a * 10 m in case 

3 

of a stress wave with (moderate) speed v ■ 10 m/sec. It is seen that E 
of quartz is of the order of, E of BaTlO^ is below, and E of Rochelle salt 

4) 

is above, the dielectric breakdown field — . While the generated voltages 

2 

V - Ea are large, the generated currents I over an area A « 1/m are only 
moderate, but the energy pulses AE are significant. Furthermore, these 
examples show that a stress polarization generator can be operated both 
below and in the breakdown regime. 

The electric polarization of solids by stress shock waves is also of 
considerable basic physical interest. Attempts at a shock wave interaction 

theory with stress polarizable solids (short circuited) were made by Alli- 
7) 8) 9) 

son — , Zeldovich — 7 , and Ivanov — 7 under various phenomenological assump¬ 
tions and compared with experimental data from piezoelectric and ferroelec- 
4) 

trie pressure transducers — . The applicability and limitations of these 
theories have been the subject of extensive discussions — — *■— , which in¬ 
dicate that the problem is not completely understood from first principles. 

Herein, a theory of the interaction of stress shock waves with (non- 
ferroelectrlc) stress polarizable solids is presented based on Maxwell's 
equations, when the external circuit contains an ohmic load. It is shown 
that the electric current pulse I(t) is determined by an inhomogeneous Vol- 
terra integral equation of convolution type. Closed-form solutions are ob¬ 
tained for the generated current I(t) by taking into consideration the dif¬ 
ferent time scales at which electric space charge and electric polarization 
by stress develop in the solid. 


40 







The following shock wave interaction theory with solids is in general 
not applicable to ferroelectric solids (except for rough estimates). Ferro- 
electrics exhibit large permanent electric polarization and a pronounced 
hysteresis effect so that the electric displacement D * D(E) is no longer a 
linear but a strongly nonlinear function of the electric field E (a hyster¬ 
esis theory for ferroelectrics is not available). 







polarization 






POLARIZATION AND CHARGE BUILD-UP 


In unstressed, shock polarizable (non-ferroelectric) solids, the electric 
dipole moments of the molecules add up to a zero vector, whereas under stress 
the dipoles are displaced preferentially in one direction so that a net polar¬ 
ization results. This is illustrated in Fig. la (o = 0) and Fig. lb (o f 0), 

which depicts the unstressed and stressed dipole configurations in a quartz 
crystal subject to a stress component o which is perpendicular to the optical 
axis. In Fig. lb, ? points to the left, but £ would point to the right if 
the crystal were rotated by 180*. For this reason, the polarization vector 
^ may point in or opposite to the direction of propagation of the stress wave 
depending on the chosen orientation of the crystal. For a stress wave o(x,t) 
propagating in the positive x-direction, the designation "±" polarization is 
used if P points into the ±x-directions. The resulting direction of the cur¬ 
rent I(t) in the external circuit is shown in Fig. 2a for the "+" polariza¬ 

tion, and in Fig. 2b for the polarization, based on Kirchhoff's law 

j Ij.dr ■ 0. 

The stress polarization generator model (Fig. 3) consists of a shock 
polarizable solid of cross section A and thickness a between plane electrodes 
at x * 0 and x * a, which are connected through an external circuit with load 
R > 0. At time t ■ 0, an explosion produced, step shaped stress wave o(x,t) 
with amplitude a and propagation speed v impinges normal onto the surface 
x ■ 0 of the solid so that, at time t > 0, the shock front is at x ■ vt and 
the stress distribution in 0 < x < a is 


o(x,t) » oH(vt - x) 


t > 0 


(1) 













where 


H(vt - x) - 1 
0 


0 < x < vt 
vt < x < a 


( 2 ) 


7 8 9) 

is the Heaviside function (compression of the slab Ax ■ a neglected) - * J 
Behind the shock front x - vt, the polarization P(x,t) builds up according 
to a relaxation process with characteristic time x. 



with 


P(x,t - 0) - 0 


(3) 

(4) 


as initial condition for the "±" polarizations (Fig. 2a,b), since the solid is 
initially unpolarlzed and unstressed. In view of Eq. (1), the solution of 
Eqs. (3) - (4) is 


P(x,t) - ±fo[l - e" <Vt ~ x)/vT ]H(vt - x) . (5) 

Accordingly, P(x,t) drops exponentially to zero for x -*• vt (shock front), and 

P(x,t) * ±f a in0<x«vt for t >> t. The polarization relaxation 

-6 4 

tine is of the order-of-magnitude x - 10 sec.— 

Due to shock polarization P(x,t) in the stressed region, which is absent 
in the unstressed region, the electric displacements behind and ahead of 
the shock front x • vt are —^ 


e l E l + p * 

0 < x < vt 

» 

(6) 

*2®;2 ’ 

vt < x < a 

9 

(7) 
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by Eqs. (5), (6), and (11), respectively. The upper and lover signs in Eq. 
(12) refer to the "±" polarization caseB. Accordingly, the space charge 
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build-up in the stressed region is determined by 


^ e - (vt " x)/vT H(vt - x) , 0 < x < vt , (14) 

31 T T T VT 

o o 


P x (x,t » 0) * o 


. (15) 


P 1 (X ’ C) * + v(t [e~ <Vt_x)/vT - e <Vt X>/VT °]H(vt - x) , (16) 


E 1 (x,t) - E^(vt,t) - (vT Q /£ l )p 1 (x,t) , (17) 

for the "±" polarization cases. The space charge relaxation time is defined 

by 

t - e /a . (18) 

o 1 

Substitution of Eq. (17) into Eq. (9) shows that the x-integration "constant" 
E^(vt,t) is determined by 


dE^(vt,t) E 1 (vt,t) I(t) 
_ + - - 

o 1 


(19) 


E x (0,0) - 0 • (20) 

Hence 

, -t/r t +t'/x 

E.(vt,t) - (Ae.) e / I(t’)e °dt’ . (21) 

1 0 


Thus, we obtain the electric fields behind and ahead of the shock front as in¬ 
tegral functionals of I(t): 


E l<x,t) 


t +t'/t 

/ I(t’)e * 
0 


•i,.» - o f 

+ e l (T_T o ) 


r -(vt— x)vt ^ 
le - e 


-(vt-xi^r 
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0 < x < vt 





vt < x < a 


(23) 


, t 

E,(x,t) * (Ae ) / I(t')dt' , 

Z Z 0 

for the "+" polarization cases. Since = 0 in the unstressed region, 

Eq. (23) obtains by direct integration of Eq. (10). 

The surface charge density p*(vt,t) » n* [D]^^ C ^ e s ^ oc ^ ^ ront is by 
Eqs. (5), (6), (7), (22), and (23) 


p (vt,t) 


(Ae,)* 1 / I(t’)dt’ - (Ae.) 
* 0 x 


-1 


-t/T 


7 I(t')e 
0 


+t'/t 


J dt’ 


(24) 


It is obvious that in general p*(vt,t) # 0, i.e., p*(vt,t) ■ 0 only if 

and T q 0. Therefore, the assumption D^(x,t) ■ D2(x,t) for x * vt made in the 

theory of piezoelectric and ferroelectric pressure transducers is in general 

incorrect. 
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INTEGRAL EQUATION FOR I(t) 


The electric current pulse I(t) in the external circuit is determined 
by Kirchhoff’s law £ E*dr * 0 for the closed current loop with load R in 
Fig. 3, where E(x,t) in the solid is given by Eqr. (22) and (23) and /t R -dr * 
+ IR. Thus, we obtain for the transient current I(t) the inhomogeneous 
Volterra integral equation: 


vt 

Ae 


t -(t-t’)/T o 

J e 

1 0 


I(t')dt’ + 


O f VT 

o 



-J— / KOdf + RI(t) 
^2 0 


-t/r 
re - re 


-t/r 


<1 - 


(25) 


The upper/lower signs refer to the "±" polarization cases. A more general re¬ 
presentation of this integral equation is achieved by introducing the dimension¬ 
less variables 

J(i) - I(t)/I 0 , l o = atff/e^A ; t - t/t Q , t Q » t ; (26) 

and 

a - a/vr , R ■ R/(o ^vt/a) , e ■ e l^ e 2 ’ 6 " T </ T ’ ( 27 ) 

as characteristic dimensionless parameters (c 1 vt/a - internal resistance of 
solid slab of thickness vt traversed by shock during polarization relaxation 
time t). By means of Eqs. (26) and (27), the determination of J(t) is reduced 
to the initial-value problem: 

r J\j(*')d V + Rj(t) 

0 0 

- ±6 (1 - ~ i~Jf -) , ( 28 ) 
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j(-t»o) - o 


(29) 


By means of the Laplace transformation and the convolution theorem, the integral 
equation (28) is reducible to an inhomogeneous differential equation of first 
order for J(s) ■ L[J(-t)] with variable coefficients of not more than third order 
in s. The general solution can then be obtained by inversion, I(t) * L ^[J(s)]. 

It should be noted that in the derivation of the integral equations (25) 
and (28), the current I(t) has been assumed to flow in the +x - direction (Fig- 3). 
Accordingly, one exspects solutions I(t) < 0 for "+"polarization and solutions 
I(t) > 0 for polarization. 
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ANALYTICAL SOLUTION 


For electrical power generation, the main interest is in shocked solids 
with sufficient pressure ionization so that 6 = x q /x << 1. For this reason, 
analytical solutions are derived for the case of small space charge relaxation 
times compared with the polarization relaxation time x. The solution for 
large space charge relaxation times, 6 = x q /x >> 1 , can also be given in closed 
form m 

For fast space charge relaxation, x q << x, the initial-value problem (28) 
- (29) reduces for the polarizations P in (+) and opposite (-) to the direction 
of shock propagation to (6 << 1): 


c~~± / + (R + t) Jtf) - ±6(1 - e~) 

5 0 


, (30) 


J(-t-O) - 0 


. (31) 


Eq. (30) consists of a cumulative voltage term, the ohmic voltage loss, and a 
delayed stress polarization voltage source. 

By means of the transformation, 


J(t) - [*(*)« 


(e/6) *] , *(t-0) 


the integral equation problem (30) - (31) is reduced to an initial-value problem 
for a first order differential equation: 











(R + i)|£ + f (R + a)<l> = ±6(1 - e -;t )e' (e/6);t 


, (33) 


ifi (-t=0) = 0 


. (34) 


Eq. (34) implies dt(/(£=0)d£ ■ 0 by Eq. (33). The solutions of Eqs. (33) - (34) 
are for the "±" polarization cases 


-f(*+R) t f(a+R)-l y , ,..y 
</<(-£) = ±6(R + t) f (R + t) (1 - e )e" t / ; cit' , 0 < t < a. 

0 

(35) 

The corresponding current solutions J(-t) for the "±" polarization cases are by 


Eq. (32): 


-*- - --- , 

5R .f(a+R)+l . 

(1+f) 6 e- (e/6);C 


0 < t < a . (36) 


These solutions can be extended to times t > a. beyond the transit time t * a if 
the assumption is made that the shock wave is not reflected at the electrode x * a. 

Eq. (36) indicates that the integral term decreases much faster with in¬ 
creasing resistance R than the ordinary stress polarization term. The 
polarization current is - by magnitude - as large as the "-"polarization current, 

(J(?)| - |J($)| , 0 < t < a . (37) 

Since the integral term is dominant by magnitude, the solutions of Eq. (36) are 
J(-t) 5» 0 for the ± polarizations and times t < a. 
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At the transit time t - a., the integral term in Eq. (36) vanishes, and 
at this particular instant the current is given by 

J (t~a) = ±6(1 - e -a )/(R + a) . (38) 

The solution in Eq. (36) shows that the (dimensional analysis) estimate 

—1 2 

I - APv/a (Table I) is correct for t < a/v and R * o a/A ~ R, since I ~ (ea/R )I 
2 

* (e^/e 2 )( R a / R ) APv/a by Eq. (36). For power applications external resistances 

R - R of the order of the slab resistance are of main interest, 
a 

Fig. 4 shows computer calculations of current pulses ^J(-t) ("±" polariza- 

_2 

tions) versus time 0 < t < a. for R a 1 and CL « 1,2, ... 10 (6/e ■ 10 ) based on 

Eq. (36). The pulse height |J( t ) j ma-r increases considerably with increasing 
a * a/vt, i.e., if the slab width & is made larger and larger compared with the 
distance vr travelled by the shock wave during the relaxation time t (saturation 
for a >> vt). 

The generated currents J(-t) are large for R $ 1, but decrease considerably 
(about R S with increasing R > 1. The steep rise of the current J(-t) within a 
time t s 10 , i.e., t <, t/ 10 is remarkable. Apparently, the complete build-up 

of the electric polarization field is not required for the development of the 
maximum current flow J(-C) if t << x o , since the electric polarization fields are 
already very large for t ~ t/10. 












CONCLUSION 


A one-dimensional theory of a solid state power generator has been pre¬ 
sented, in which the electromotoric force is produced by electric polarization 
of a solid in the stress field of a shock wave. The numerical results indicate 
that energies of the order AE - 10^ Joule can possibly be generated per pulse. 
This type of power generator does not require an external magnetic field and 
appears to be promising as a pulsed power source for special applications 
(for which the efficiency of transformation of chemical energy of explosives 
into stress waves and electrical energy is of no concern). 

Experiments should be carried through to verify the theoretical expec¬ 
tations. In comparing experimental and theoretical current pulse forms, it 
should be noted that the temporal current distribution depends noticeably on 
the shape of the stress wave. Current pulse forms for other stress wave shapes 
than the simple step wave can be calculated in complete analogy. It might also 
be of Interest to refine the theory by including multi-dimensional effects and 
deformation of the solid. 
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INITIAL-BOUNDARY-VALUE PROBLEM FOR DIFFUSION OF MAGNETIC FIELDS INTO CONDUCTORS 


WITH EXTERNAL ELECTROMAGNETIC TRANSIENTS 

H.E. Wilhelm 

Michelson Laboratory, Naval Weapons Center, China Lake CA 93555 


ABSTRACT 


The initial-boundary-value problem for the diffusion of an in¬ 
itially homogeneous magnetic field into a slab of conductivity o < 00 
and width Ax * 2a is solved, under consideration of the electromagne¬ 
tic wave pulses generated at the surfaces of the conductor by its 
interaction with the external magnetic field, which propagate into 
the surrounding vacuum. The analytical solutions show that (i) the 
external electromagnetic transients are necessary in order to correc¬ 
tly satisfy the boundary conditions for the tangential electric and 
magnetic field components, and (ii) the spatial and temporal develop¬ 
ment of the electromagnetic field and electric current in the conduc¬ 
tor is quantitatively determined by a new dimensionless parameter group 
- 1/2 

R “ U Q aac [c ■ (^ 0 e 0 ) 1 * This "magnetic Reynolds number of the va¬ 

cuum" determines the coupling between the transient fields in the con¬ 
ductor a > 0 and the ambient space (a ■ 0). 


57 





INTRODUCTION 


The diffusion of electromagnetic fields B(r,t), E(r,t) in a conductor of finite 

conductivity a and normal surface vector n(s), when the electromagnetic field ^ Q (r,t) 

and ^ (r,t) outside of the conductor are known, is in general described by Maxwell's 
o 

equations without displacement current, where the tangential field components are 

assumed to satisfy the boundary conditions^’ 2 ^ n *[B(r,t) - B Q (r,t)] * t) and S * 

f£(r,t) - E o (r,t)] • If the external electromagnetic field is time-independent 

and electric potential fields are absent, then B q * B Q (r) and E q * 0 (since V x E q * 

-3B /St * t) and E * -V$ 2 0), so that the tangential boundary conditions reduce to 1- ’ 2 ^ 

o o o 

n x [B(r,t) - B Q (r) ] * <5 and n x E(r,t) * T,. These boundary conditions have found wide- 

1 ) 2 ) 

spread use in mathematical physics , eleccromagnt ic theory , and the theory of mag- 

3 4) 

netic flux compression (at the outside surface of the liners) ’ . However, these 

boundary conditions are questionable approximations, since they do not take into con- 

slderation the wave fields B(r,t), E(r,t) propagating away from the conductor into the 

-1 

surrounding medium, which have their sources in the transient current fields j * 

5 k B of the conductor. 

For a concrete illustration of the problematics, consider the diffusion of an 
external (homogeneous) magnetic field, B q ■ {0,B o> 0} for )x| > a, into a conducting 
slab in the region |x| < a which is field free at time t * 0 (Fig. 1). Using the con¬ 
ventional boundary conditions, the transient magnetic field B(x,t) * {0,B(x,t),0) in 
the conductor is determined by the parabolic initial-boundary- value problem:^ 

3B/3t - k3 2 B/3x 2 , |x| < a , t > 0 , (1) 

B(x,t-0) - 0 , |x| < a , (2) 

B(x«±a,t) * B , t > 0 , (3) 

o 


where k * l/vi Q a, By means of Fourier's method, the general solution of Eqs. (1) - (3) 
is obtained as a superposition of eigenfunctions^: 



t > 0, (4) 








Since V x B - p oE, the electric field 


with B(x,t) -*• B q in jxj < a for t -*■ °°. 
E(x,t) * (0,0,E(x,t)} in the conductor is 


2B oo ')')') 

. _o. r . _.n-l -ic(2n-l) ir t/4a . 2n-l 

E(x,t) = —— I (-1) e sin —— ttx 

^ n=l 


2a 


< a, t> 0. 


( 5 ) 


In accordance with the boundary conditions (3), the space outside of the conduc¬ 
tor remains unperturbed while the electromagnetic field diffuses into the conductor, 

B Q (x,t) « B q , E o (x,t) - 0 , |x| > a , t > 0 (6) 

The transient currents j * V x b/p q in the conductor are "eddy currents," and, 
therefore, cannot produce transient magnetic fields B Q (x,t) - B Q (x,t) - B q i 4 0 
in the outside region |x| > a. The net current I(t) through any cross section 
z ■ constant vanishes, due to the boundary conditions (3): 

g 

I(t)/Ay * p” 1 f+ a [3B(x,t)/3x]dx ■ p" 1 / ° dB * 0 (7) 

-a B 

o 


By comparing the conductor solutions (4) and (5) with the vacuum solutions (6), 

it is seen that B(x ■ ±a,t) - B * 0, but E(x * ±a,t) - E (x «* ±a,t) * E(x * ±a,t) f 4 0! 

o o 

1 , 2 ) 

Thus, the conventional boundary conditions lead to a violation of the fundamental 
law of the continuity of the tangential electric fields at Interfaces. 

The correct formulation of the boundary conditions requires consideration of 
the simultaneous wave fields B + (x,t), E + (x,t) propagating with the speed of light c 
into the positive and negative half-spaces x > +a and x < -a surrounding the 
conductor (Fig. 2), which are excited by the transient current fields j(x,t) - 
P Q ^3B(x,t)/3x in the (space-charge free) conductor. No matter how small the external 
transients B ± (x,t) and E ± (x,t) are (in comparison with B q T 4 0 and E ■ 0), they have 
to be taken into account in order to rigorously satisfy the boundary conditions 
n x [$] •• 5 and n * [E] ■ 5 for the continuity of the tangential electromagnetic 
fields at conductor interfaces. 
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The quantitative assessment of the significance of the external transients of 
the diffusion process leads to the discovery of a new dimensionless parameter combi¬ 
nation, which has the physical meaning of a "magnetic Reynolds number of the vacuum": 

R * ii Q aac , c = (i i o e o )“ 1/2 * 3 * 10** m/sec . (8) 

In the following, the formulation of the initial-boundary-value problem for 
diffusion processes with external transients and its analytical solutions for the 
transient electromagnetic fields inside and outside the conductor are presented. The 
qualitative and quantitative importance of the new boundary conditions and the 
external wave fields are discussed in terms of R. 

The presented theory has important implications for the evaluation of the flux 

3 4) 

losses through the liners of magnetic field compressors ’ , the electromagnetic 

acceleration of conducting macroparticles** ’ the electromagnetic induction in 

8 9) 

conductors moving relative to external magnetic fields ’ ', and for the interaction 
of transient plasma shock waves with external magnetic fields 10,11 ^. The general 
significance for theoretical physics is obvious. 
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INITIAL-BOUNDARY-VALUE PROBLEM 


Subject of the consideration is the diffusion of the magnetic field into a con¬ 
ducting slab |x) < a, which is initially embedded in a homogeneous magnetic field 

B * (0,B ,0,}, under sinultaneous emission of electromagnetic waves from the conduc- 
o o 

tor surfaces x * ±a (Fig. 1). 

The transient electromagnetic fields 1t + = {0,B + (x,t),0} and * (0,0,E + (x,t)} 
in the infinite vacuum half-spaces (o » 0, y a y Q ) x > +a and x < -a are determined 
by the initial-boundary-value problems (±) for the wave equation. 


aV/3t 2 


2 2 2 
cl B./ax^, 


±x > a. 


t > 0, 


(9) 


B.(x,t * 0) * B , ±x > a , (10) 

± o 

B + (x » ±a,t) » B q + ij» + (t) , t > 0, (11) 

since 

3E + /31 * c 2 3B + /3x , 3E + /3x = 3B + /3t (12) 

by Maxwell's equations with displacement current. The solutions of Eqs. (9) - (11) 
for the still undetermined boundary values i{i + (t) are: 


and 


B ± (x,t) - B q + * ± (t +~ ) 


- B 


o 


a < ±x < a + ct 
a + ct < ±x < » 


E + (x,t) - +c<j/ + (t + ) 

- 0 


a < ±x < a +■ ct 
a + ct < ±x < • 


(13) 


(14) 


These solutions are typical for hyperbolic equations, i.e., the boundary values i|< + (t) 
are "transported" into the half-spaces ±x > a with the speed of light c, so that dis¬ 
continuous wave fronts result at x ■ ±(a + ct). 
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Let the external magnetic field B q be switched on at t = 0 so that no electro- 

1 2 ) 

magnetic fields exist in the conductor for t < 0 * . The conductor has a finite con- 

12 ) -*■ ■+• 

ductivity a and can, therefore, not carry surface current densities, j* ■ ^ij^oEAx = 0 
for o < 00 and E bounded. Accordingly, the boundary conditions for the tangential electric 
and magnetic field components at the conductor vacuum interfaces are 


B(x » ±a,t) - B q + i|/ ± (t) 
E(x - ±a,t) = +ci|» ± (t) 

where 

E(x,t) * ic3B(x,t)/3x 

by Ohm's law is the electric field 
in |x| < a. Furthermore 


t > 0 , 


(15) 

t > 0 , 


(16) 

M < a 

t > 0 , 

(17) 


in the conductor, and B(x,t) is the magnetic field 


k * 1/u a > 0 . (18) 

o 

(The boundary conditions n-[eE] = 0 and n*[B] * t) are satisfied since ^ and 1$ have no 
normal components.) By elimination of the unknown boundary values E(x * ±a,t) and <|i + (t) 
from Eqs. (15) - (17), boundary conditions involving only the magnetic field B(x,t) in 
the conductor are obtained: 


(19) 


These are the new boundary conditions for the diffusion of magnetic fields B(x,t) into 

1 2 ) 

conductors. They differ from the conventional boundary conditions * through the curl 
terms 3B(x ■ ±a,t)/3x 4 0, which consider the emission of magnetic dilution waves from 
the conductor surfaces x ■ ±a into the vacuum | x| > a. 

Within the conducting slab of finite width 2a, the propagation of the magnetic 

1 2 . 12 ) 

field can be treated in the nonrelativistic or diffusion approximation ' . Accor¬ 

dingly, B(x,t) in the initially field free conductor is determined by the parabolic 
initial-boundary-value problem: 
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2 2 

3B/3t * «3 B/3 x , I xl < a , t > 0 , 


B(x,t * 0) = 0 


3B(x = ±a,t)/9x ± hB(x=±a,t) = ihB^, t > 0 


h * c/k > 0 


The transformation. 


B(x,t) * B + B(x,t) , |x| < a , t > 0 , 


reduces Eqs. (21) - (22) to an initial-boundary-value problem with standard "radiation" 
boundary conditions: 

3B/3t * k3^B/3x^ , |x| < a , t > 0 , (25) 

B(x,t - 0) = -B q , |x) < a , (26) 

3S(x = ±a,t)/3x ± hB(x=±a,t) =* 0 • (27) 

In accordance with Fourier's theorem, the solution of Eqs. (25) - (27) is obtained 
as a superposition of eigen-functions B n (x,t) for the region |x| < a which satisfy the 
boundary conditions (27): 

00 (h^+k^)k ^sink a -xk^t 

B(x,t) - -2B l -S—2-— e n cosk x , |x| < a , t > 0 , (28) 

n-1 [ (l/-H^)a-Hi] n 

where L 

k n a tg(k Q a) - ha , n - 1,2,3,..., (29) 

gives the eigenvalues ^associated with the boundary conditions (27). 

1. Conductor Solutions: 


For a compact representation of the field solutions, dimensionless independent 


(£,t) and dependent variables are introduced, 
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5 = x/a 


(30) 


2 

t = Kt/a 


a = k a 
n n 


8(5,t) = B(x,t)/B o , Eft, t) = E(x,t)/(KB o /a), /ft,x) = j(x,t)/(B o /p o a). 

(31) 

According to Eqs.(24) and (28), the solutions for the dimensionless fields 8ft, t). 
Eft, t) =* 38ft,x)/3£, and /ft, t) in the conductor are given by: 


,-02 2 . -1 . 2 
• (K +a )a sina -a t 
r _n n_n n 


Bft ,x) 

„ „ v nn n n 

“1-22, ,9 e coso E . 

n=l [(R z +x/)+R] n 

Id < 1 , t> 0 , 

(32) 

E(C.t) 

/»2. 2. . 2 
® (R +o ) sim -a x 

= —7^1—“ e n sin v 

n«l [(R +a )+R] 
n 

|5| < 1 , t > 0 , 

(33) 

/ft. t) 

” 8(5,t) , 15! 1 1 . r >_ 0 , 


(34) 

» tg& 
n n 

» R , n « 1,2,3,..., R * ac/< * p Q oac 


(35) 


_2 

by Eqs. (8) and (29). For sufficiently large times x >> , the homogeneous magnetic 

field has diffused completely into the conductor. 


8ft, t) -*• 1 , Eft.x) -*■ 0 , /ft,x) -*■ 0 , t -+• 


(36) 


In the hypothetical limit of infinite magnetic vacuum Reynolds number R, Eq. (32) re¬ 
duces to the known solution (4) for the conventional boundary conditions'’), 


. 4? (-I) 11-1 -(2n-l) 2 ir 2 T/4 2n-l . , _, „ . 

limBft ,t) - 1 - — l e cos— 7~ |5| < 1, t > 0, 

R-*» n«l v ” ' 


since 


lim o ■ 2 °^^~ tt, n ■ 1,2,3,... 
R-*» n 2 


(37) 

(38) 


Comparison of Eq. (32) with Eq. (37) indicates that the 8ft,x) solutions 


with the new and conventional boundary conditions differ not much if R >> * it/2. 
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2. Vacuum Solutions: 


In view of the boundary conditions (15), Eq. (32) yields for the boundary 
values ¥ (t) = 6(£=±1,t) - 1. Accordingly, Eqs. (13) and (14) give for the electro¬ 
magnetic fields in the positive (£ > +1) and negative (£ < -1) half-spaces the so¬ 
lutions: 


B + (5,t) - 1 + ¥ ± (t + ^) , 1 < ±£ < 1 + Rt , 


= 1 


, 1 + Rt < ±i < 


and 


£ ± (5,x) = *R¥(t + ^) 


= 0 


1 < ±£ < 1 + Rt , 
1 + Rt < ±g < ® , 


where 


Mt = - 


(R^+n^)sina n cosa n -a^[r+(£+l)/R] 


n=l a. [ (R 2 +a 2 )+R] 


(39) 


(40) 


(41) 


Eqs. (39) and (40) represent electromagnetic wave pulses which propagate with the 
dimensionless spaed R(c) from the conductor surfaces £ • ±1 into the vacuum spaces 
> 1 with discontinuous wave fronts at £ = ±(1 + Rt). They are kicked on at t = 0 
and their emission lasts to the end (t -► ®) of the diffusion process in the conductor. 
The vacuum fields 8 + (?,t) are in opposite direction of B q , i.e., they represent di¬ 
lution waves (Fig. 2). 

In the case of large coupling numbers, R >> 1, Eqs. (39) and (40) yield with 


Eq. (41) 


8 ± (C,t) * 1 + 0[R _1 ] 
E + (C.x) *±2 l e 2 

n-1 


,2n-l ,2 
-(' r— it) t 


1 < ±£ < 1 + Rt , t > 0, R » 1, (42) 
1 < ±£ < 1 + Rt , t > 0, R >> 1, (43) 


since cosa * (-l) n *(2n-l)W2R for R >>mrby Eq. (35). 
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The magnetic field B q outside of the conductor remains nearly unperturbed 
during the diffusion process, B + ~ R \ whereas the external electric transients 
4 0 are of order R^ behind the wave fronts, f, = ±(1 + Rt), for R >> 1. However, 
since V*B + = c ^3E + /3t, not only - R but also 8 + - R ^ cannot be neglected for 
R >> 1. 

Thus, we have shown how self-consistent solutions can be obtained for the 
electric and magnetic fields in the conductor and the surrounding vacuum, which 
satisfy the boundary conditions for the continuity of the tangential electric and 

magnetic fields at the conductor-vacuum interfaces. The conventional boundary 

1 2 ) 

conditions for electromagnetic diffusion processes * , ignore the external elec¬ 

tromagnetic transients, violate the boundary condition for the tangential electric 

field, and permit no Poynting vector S = ExH outside of the conductor. As a result, 

12 ) 

the conventional boundary conditions make it impossible for electromagnetic 

fields to diffuse through condi:~c f and to escape into the ambient space. 

For both the conductor tc ..':m solutions, the limit R ■+■ 0, which implies 
a -*■ 0 since a 4 0, is not realizable since the conductivity of conductors is by 

definition large. For insulators or extremely poor conductors (R * a -*■ 0), the 

, , 12 ) 
nonrelativistic or parabolic diffusion equation is inapplicable . Therefore, 

the investigation of the limit R -*• 0 would require solution of Maxwell's equations 

with displacement current in the slab |x| < a given elsewhere^^. 

The generalization of the theory required for conductors and external media 

(vacuum, gases, fluids) with different permittivities e and u are trivial but 

complicate the notation. 
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DISCUSSION 


It is known that Maxwell's equations with displacement current and the non- 
relativistic Ohm's law j = oE combine to a hyperbolic diffusion equation for the 
magnetic field B(r,t) in conductors'^ , 


3 2 ^ , 1 3 $ 2 „ 2 ± 

where 


(44) 


T R - e Q /a (45) 

is the field relaxation time, e.g., x 0 a (10 ^/36 it)/ 6*1Q 2 * 10 ^sec for copper 

K 

with 0 = 6*10^0 ^/m. Eq. (44) reduces to the parabolic diffusion equation in the 
limit t r << 1: 


3B 

3t 


2_2 1 
z V B 


« 1 


(46) 


The parabolic diffusion equation is an excellent approximation, since the relaxa¬ 
tion time of conductors is extremely small, x^ << 1, By Eqs. (45) and (46), the 

field relaxation time x and the diffusion time x are interrelated by 
K D 


T 


-l 

D 


2 , 2 

: /a 


2 

r * Mia 
D o 


(47) 


where a is the extension of the conductor. For conductors, the diffusion time is 

-7 7 -4 

relatively large if a^ is not microscopically small, e.g., x^ * 4**10 6*10 *10 

-2 -2 
~ 10 sec for a copper slab of width a * 10 m. 

2-+ 2 

Comparison of the neglected term 3 B/3t with the second and third terms of 

Eq. (44) reveals the relation of the parabolic diffusion approximation to the new 

coupling number R * p oac: 

o 


6b 




-i alj 
R 3t ! 


at 


: 2 V 2 b! 


1r . Vo „ r -2 

T D (u aa) 2 
o 


(48) 


This result again confirms the validity of the parabolic diffusion equation for 

conductors, for which R = u oac >> 1. E.g., R = 4itx10~ 7 6x10 7 x10~ 2 3x10 8 - 10 8 

o 

_2 

for a copper slab a = 10 m. More important, Eq. (4g) demonstrates that the 

2 -*- 2 -2 
neglected relativistic term 3 B/3t in the conductor is small of order R <<<< 1, 

whereas the calculated electromagnetic fields in the conductor are of order 

8 - E - R® [Eqs. (32) - (33)], and the external electromagnetic transients are 

of order 8 + - R ^ and E + - R^ [Eqs. (39) - (40)], since in Eq. (41) for large R 

| cosaj - [1 +(R/a^) 2 ] ^ 2 * a fl /R - -^2——^rrR 1 , R » nir , n » 1,2,3,... . 

(49) 

In conclusion, it is noted that, in conductors, magnetic field diffusion is 
a nonrelativistic process (as is electric conduction, j * aE). The electric tran¬ 
sients E + in vacuum must be of the same order as the electric field E in the con¬ 
ductors, E + ' E ~ R^, since otherwise the tangential electric field could not be 
continuous across the conductor-vacuum interface. On the other hand, the external 
magnetic transients 8 + are small of order R 1 * (y Q oa) ^c 1 since the magnetic 
field energy flows with the speed of light towards the conducting cavity. The 
deeper physical reason for these electromagnetic transients is to be seen in the 

conservation laws for electromagnetic energy and momentum, which follow from Max- 
12 ) 

well’s equations . 
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INITIAL-BOUNDARY-VALUE PROBLEM FOR ELECTROMAGNETIC INDUCTION IN ACCELERATED CONDUCTORS 

MOVING ACROSS MAGNETIC FIELDS 


H. E. Wilhelm 
Michelson Laboratory 
Naval Weapons Center 
China Lake, California 93555 


Abstract 

Boundary conditions are derived for the Interfaces of a conductor moving across 
an external magnetic field in an ambient medium (vacuum or non-conductor), which 
consider the emission of electromagnetic waves from the conductor surface as a 
result of electromagnetic induction. These boundary conditions are applied to 
the initial-boundary-value problem for the electromagnetic induction in a con¬ 
ducting slab, which is accelerated across a homogeneous magnetic field to a non- 
relativistic velocity. Fourier series solutions are presented for the transient 
electromagnetic fields in the moving conductor and the discontinuous electromag¬ 
netic waves in the ambient space. It is shown that the transient electromagnetic 
fields inside and outside the conductor are due to two mechanisms, i.e., "velocity 
induction" (ordinary induction) and "acceleration induction" (dv(t)/dt/O). 
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INTRODUCTION 


Although the theoretical foundations for electromagnetic induction in con¬ 
ductors moving across magnetic fields were formulated in 1908 by H. Minkowski^, 

only simple problems such as stationary unipolar induction in rotating discs 
2 ) 

have been discussed. The electrodynamics of moving media has been a subject 

3) 

of basic research and also controversy to date. In the treatment of electro- 

4 5) 

magnetic induction, e.g., in liners of magnetic field compressors ’ and tran¬ 
sient plasma shock waves interacting with external magnetic fields B it 

has become customary to use the boundary condition for the tangential magnetic 
field in the form nx[ B - B q ] * ? at the conductor-gas interface, where B is 
the transient magnetic field in the conductor and B q is the unperturbed (!) 
external magnetic field. This boundary condition leads to electromagnetic 
solutions B,E ir. the conductor which do not satisfy the corresponding bound¬ 
ary condition ux{E - ? q ] = 5 for the tangential electric fields at the mov¬ 
ing interface. These "conventional" boundary conditions produce approximate 
to incorrect results, depending on the physical situation. 

For an illustration of the problematics of the conventional boundary 
4 5) 

conditions, ’ which are also being used in the analysis of magnetic field 

8 ) 

diffusion into conductors at rest , consider a conducting slab Ax * 2a with 
its surfaces initially at x ■ ±a in a transverse homogeneous magnetic field 
- {0 ,B q , 0} for |x| <» (Fig. 1). At time t ■ 0, this conductor is set 
into motion with a velocity v(t) ■ {di(t)/dt,0,0} so that its front and rear 
surfaces are at x * x(t) i a at time t > 0 where £(t*0) * 0. No matter whether 
the induction of the transient magnetic field iT(x,t) in the moving conductor 
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9) 

is described by the relativistic wave equation or the non-relativistic dif- 
9) -f ^ 

fusion equation , the initial condition B(x,t=0) = B q , -a < x < +a, and the 
conventional boundary conditions B(x=x(t)±a,t) = B q , t > 0, permit only one 
and the same solution, B(x,t) * B q , which implies E(x,t) * -v(t)xB Q , and 
J(x,t) ■ 5 by Ohm's law for the moving conductor. These simple solutions are 
due to the conventional boundary conditions without external perturbations and 
are obviously not correct, since the boundary condition for the tangential elec¬ 
tric field is not satisfied, nx[i£ - E ] * v(t)B e # 5 where it * Z. Since M. 

o o y o 

Faraday it is an experimental fact that transient electromagnetic fields and 
currents are induced in the conducting slab as soon as it is moved relative to 
the external magnetic field B q . 

In the following, an analytical solution is presented for the initial- 

boundary-value problem of the electromagnetic induction in a conducting slab 

Ax = 2a, which is at rest for t < 0 in a transverse, homogeneous magnetic 

field B , and which is set in motion at t - 0 with a nonrelativistic velocity 
o 

v(t) of arbitrary (finite) acceleration dv(t)/dt (Figs. 1,2). The electro¬ 
magnetic induction in the conductor is shown to produce transient electro¬ 
magnetic fields within the conductor and electromagnetic waves at the moving 
conductor surfaces x “ i(t) ± a, which propagate with the speed of light in 
the surrounding space to infinity. Only if the electromagnetic waves outside 
of the moving conductor are taken into consideration, consistent solutions of 
Maxwell's equations exist which satisfy the boundary conditions for the tan¬ 
gential and normal electromagnetic fields at the conductor-non-conductor 
Interfaces. 

The analytical solutions for the moving conductor and surrounding spaces 

permit significant conclusions. The transient electric field induced in the 
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conductor is the sum of a field which is proportional to the velocity v(t) 
and a field which is an integral functional of the acceleration dv(t)/dt 
of the conductor. A fundamental dimensionless group (a * conductivity, 
permeability of conductor) 


R ■ Ujoac, c 


(u 2 e 2^ 


- 1/2 


is found which represents a "magnetic Reynolds number" of (i) "free space" 

(if the conductor moves in vacuum or gas with permeabilities £ o ,p q ,), or (ii) 
"non-conducting space" (if the conductor moves in a nonconducting medium with 
permeabilities e^v^). R determines the coupling of the transient electromagnetic 
fields inside and outside the conductor. The external magnetic transients 
would be negligible for R -*• °°, however, the external electric transients do 
not vanish for R -*■ <*>. Thus, the conventional boundary conditions which assume 
that the external magnetic field B q remains unperturbed outside of the conductor 
while transient electromagnetic processes occur in the latter are not correct 
within electromagnetic theory. 
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I 


,C4uv, 


BOUNDARY CONDITIONS 

For the analysis of electromagnetic Induction in moving conductors (conducti¬ 
vity a < °°, permittivities which move with a velocity field v(r,t) relative 

to the "laboratory system" in a non-conducting medium (vacuum, gas, or fluid; o * 0, 
£ 2 *^ 2 ^* t * le boundary conditions for the electromagnetic fields (conductor) 

and 1*2> (non-conductor) are required in the L-frame. Integration of Maxwell's 
equations with displacement currentover the interface 1-2 with velocity v(r,t) 
and normal vector n(r,t) (direction 1 -*• 2) yields for the tangential and normal elec¬ 
tromagnetic field components the boundary conditions in the L-system: 


nxf &2 ~ = +(n-v) [^ 2 - 1^] 

nx['S 2 /u 2 - ■ 1* ~ (n*v)[e 2 ^ 2 “ e^] 

n* [£2^2 ** E I^1^ * P 

n • [\ - \] * 0 


where 


V 


P v 


lim pAs 

Aa+0 

IpI-m- 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 


are the surface current and charge densities of the interface As -*• 0. A conductive 

surface current density does not exist at a conductor of finite conductivity, ^ * 

lim oSas ■ ^ for a < «• and 1 bounded. 

As+0 

In many cases, the permittivities of good conductors and their ambient atmos¬ 
pheres equal the free space values, ^ m e V 1 2 " With this simplification 
in notation, the boundary conditions (1) - (4) are applied to the front (+) and rear 
(-) surfaces x ■ x(t) ± a of a conducting slab (o < ») with the fields (Fig. 1) 
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B - {0,B(x,tXO} » E ■ {0,0,E(x,t)} , x(t) - a < x < x(t)+a ; 

B(x,t*0) = B q , E(x,t=0) - 0 , -a < x < +a, , (6) 

which moves with the velocity v(t) * {dx(t)/dt,0,0} in an ambient medium (a * 0) with 
the fields (Fig. 1) 

B + - {0,B + (x,t),0} , E + * {0,0,E + (x,t)}, x% x(t) ± a ; 

B + (x,t“0) * B q , ^ + (x,t-0) * , x % ±a . (7) 

B^is an external homogeneous magnetic field which fills uniformly the conductor (1) 
and the medium (2). The boundary conditions (3) and (4) are satisfied identically 
since B,E and B + ,E + are 1 n so that "j * 0 and p - 0 by Eq. (5). The tangential 
boundary conditions (1) and (2) yield for the fields (6) - (7): 

E + (x»x(t)±a,t) - E(x-x(t)±a,t) » +v(t) [B ± (x»x(t>±a,t) - B(x-=x(t)±a,t)], 

" (S) 

B + (x=x(t)±a,t) - B(x*x(t)±a,t) * + v(t)c~ 2 [E + (x=x(t)±a,t) - E(r-x(t)±a,t)], 

(9) 

where 

c - (ue)" 1/2 (10) 

is the speed of light. For non-relativistic conductor motions, Eqs. (8) - (9) reduce 
to: 

E + (x*x(t)±a,t) - E(x-x(t)±a,t) - 0 , v(t) 2 « c 2 , (11) 

B + (x-x(t)+a,t) - B(x-x(t)±a,t) - 0 , v(t) 2 « c 2 . (12) 

According to the non-relativistic Ohm's law for conductors which are moving with 
a velocity v relative to the L-system, 1 ■ o(f + vx$), the electric field ?(x,t) in 
the conductor is expressed in terms of 3(x,t), 
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E * -vB + (pc) ^3B/3x , x(t) - a < x < x(t) + a , v(t) 2 << c 2 . (13) 

The electromagnetic field 1i + (x,t), E + (x,t) in the adjacent semi-infinite half spaces 
is determined by the hyperbolic initial-boundary-value problems: 

3 2 B + /3t 2 - c 2 3 2 B + /3x 2 , x ^ x(t) + a , t > 0 , (14) 


B.(Xjt^O) * B 


, x % ±a 


B + (x«x(t)±a,t) = B(x-x(t)±a,t) , 


t > 0 , 


3E + /3x =» 3B + /3t , 3B + /3x - c 3E + /3t 


by Maxwell's equations for homogeneous non-conductors. Eq. (16) couples the solu¬ 
tions B + (x,t) in the semi-infinite spaces x * x(t) ± a to the solution B(x,t) in the 
conductor, x(t) - a < x < x(t) + a. By Eqs. (14) - (17), the ambient electromagnetic 
field transients are of the form: 


B + (x, t) » B o +? t (t + X „ a ) , x(t) + a ^ x •> ± (a + ct) , 


, a + ct < ±x < ® 


E + (x,t) - +cf + (t + ~ ■— ) , x(t) ± a ^ x ^ + (a + ct) . 


, a + ct < ±x < ® 


V± ( t + 2^-) - B(x-x(t)±a,t) - B q (20) 

by Eq. (16) determines the form of the wave functions ¥ + 0l + ) of the "self similar" 
arguments n ± * t + x/c + a/c from the conductor solution B(x,t). 

The solutions (18) - (19) are typical for hyperbolic initial-boundary-value 
problems, i.e., the boundary values V (t + x(t)/c) are transported with the speed of 


L 




light c into the half spaces x «? x(t) ± a, so that discontinuous wave fronts result 
at x ■ ±(a + ct). 

By means of Eqs. (18) - (20), the boundary-values B + (x=x(t)±a,t), E + (x=x(t)±a,t), 
and ¥ (t £ x(t)/c) are eliminated from the boundary conditions (11) and (12): 

~~ ^i f x ± * ,t) ' 1 t c + v(t)]B(x=x(t)±a,t) = ±cB q , v(t) 2 « c 2 . (21) 

This are the fundamental new boundary conditions for moving conductors which (i) in¬ 
volve only boundary values of the magnetic field B(x,t) of the conductor, and (ii) 
consider the emission of electromagnetic waves, 

S ± (x,t) - ¥ ± (t + 2~-^) , £ ± (x,t) - +cf(t + (22) 

from the conductor surfaces x *■ x(t) ± a into the ambient spaces x x(t) ± a. 

Since the magnetic field in the conductor is the sum of the external B q and a tran¬ 
sient B(x,t), 

B(x,t) * B q + B(x,t) , x(t) - a < x < x(t) + a , (23) 

Eq. (21) gives for the transient conductor field the boundary conditions: 

± [ c + v(t)]B(x"x(t)±a,t) - v(t)B Q , v(t) 2 << c 2 . (24) 

ii 4 -1 

In view of the difficulties to accelerate macroscopic bodies to speeds )v| > 10 m sec , 
the stronger non-relativistic condition jv| << c is in general fulfilled, which redu¬ 
ces Eqs. (21) and (24) to the simpler boundary conditions: 

(po)" 1 ± c[B(x*x(t)±a,t) - B o J - v(t)B o , jv(t)| « c , 

(25) 
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and 



(pg)" 1 ± cB(x=x(t)±a,t) = v(t)B Q , jv(t)| « c . (26) 

If the Vxli and v(t)B Q terms are omitted, Eqs. (25) and (26) reduce to the con- 
4-7) - - 

ventional boundary conditions ' B(x=x(t)±a,t) = B q and B(x*x(t)±a,t) = 0, respec¬ 
tively. Comparison shows that the conductor currents VxB/p and the induced currents 
av(t)B Q at the conductor surface are the sources of the emitted electromagnetic 
transients. 




INITIAL-BOUNDARY-VALUE PROBLEM 


Consider a rigid conducting slab of width = 2a with surfaces in the planes 
x = x(t) ± a at time t > 0 (Fig. 1). This conductor is exposed to an external mag¬ 
netic field = {0,B o> 0} which is homogeneous throughout the space -« < x < +», and 
is accelerated to a (non-relativistic) velocity v(t) = {v(t),0,0} from an initial po¬ 
sition at rest x(t=0) = 0 so that 


with 


v(t) * dx(t)/dt , x(t) * / v(t')dt' , dv(t)/dt £ 0 

0 


x(t*0) = 0 , v(t=0) * 0 , dv(t=0)/dt £ 0 


, (27) 

• (28) 


The electromagnetic induction of the transient electromagnetic fields $ * {0,B(x,t),0}, 
? = {0,0,E(x,t)> in the conductor of finite width, as a result of its accelerated mo¬ 
tion v(t) across the external magnetic field 5 q , is determined by the parabolic ini¬ 
tial-boundary-value problem for B(x,t): 


-~7 + v(t)~■ = , x(t) - a < x < i(t) + a , t > 0 

3C 3x 3x Z 


B(x,t-0) - B 


-a < x < +a 


(29) 

(30) 


3 . B (x*x l t) ,± a,t) ± h { B (x-x(t)±a,t) - Bl - <~ 1 v{t)B , t > 0 . (31) 

OX O O 

where 

k * 1/uo , h ■ c/k . (32) 

Eq. (29) follows from Maxwell's equations and Ohm's law j = cr(E + v*B) for moving con- 

9) 

ductors in the diffusion approximation, * e/o -► 0. The boundary conditions (31) 
couple the electromagnetic induction process in the conductor to the external transi¬ 
ents in the ambient medium [Eq. (25)]. The transformation. 


B(x,t) - B q + B(x,t) , -a < X < +a , t>0 
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( 33 ) 







X * X - x(t) 


t > 0 


( 34 ) 


reduces Eqs. (29) - (31) to the initial-boundary-value problem: 

- 2 ~ 2 

3B/31 - <d B/3x , -a < x < +a , t > 0 , (35) 

B(x.t=0) - 0 , -a < x < +a , (36) 

3B(x*±a,t)/3x i hB(x=±a,t) = K -1 v(t)B Q , t > 0 , (37) 

The linear IBVP (35) - (37) is decomposed into a BVP and an IBVP by means of the 
ansatz: 

B(x,t) * F(x,t) + G (x, t) , -a < x < +a , t > 0 , (38) 

where 

3 2 F/3x 2 * 0 , -a<x<+a , t>0 , (39) 

3F(x*±a,t)/3x ± hF(x*±a,t) - K -1 v(t)B o , t > 0 , (40) 

and 

3G/3t - k3 2 G/3x 2 - 3F/3t , -a < X < +a , t > 0 , (41) 

G(x,t»0) - -F(x,t»0) , -a < x < +a , (42) 

3G(x*±a,t)/3x ± hG(x=±a,t) « 0 , t > 0 . (43) 

The solution F(x,t) of Eqs. (29) - (30), the source 3F(x,t)/3t in Eq. (41), and the 
initial value F(x,t-0) in Eq. (42) are: 

F(x»t) - [v(t)B Q /2(K + ac)]x , -a < x < +a , t > 0 , (44) 

3F(x,t)/3t ■ I ^jt^ b 0 /2 ( k + ac )]x * -a < x < +a , t > 0 , (45) 

and 

F(x,t-0) - 0 , -a < X < +a (46) 
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since v(t”0) * 0- With 3F(x,t)/3t odd in X and F(x,t=0) = 0, the initial-boundary-va¬ 
lue problem (41) - (43) is solved by means of the Fourier expansions: 


G(x, t) = I G (t)sin k x , -a < x < +a 
", n n 

n*l 


where 


n«l 


k a ctg k a * -ha , n * 1,2,3, ... 


yields by means of the inverse Fourier theorem 


dG (t)/dt + <k^G (t) - -S (t) 
n n n n 


where 


G (t-0) - 0 


S (t) * (a/«)K B dv(t)/dt 
n no 


K * (h^ + k^)sin k a/(k a)^[(h^ + k^) + (h/a)] 
n n n n n 


G (t) - - / e 
0 


t -Kk*(t-t') 


S (t')dt' 
n 


t > 0 , 

(47) 

t > 0 , 

(48) 

* 

(49) 

t x associated with 

the 

n 


) into Eqs. (41) and 

l (42) 

t > 0 , 

(50) 


(51) 

t > 0 , 

(52) 

• 

(53) 

is 


t > 0 

(54) 


Combining of Eqs. (47), (52), and (54) yields as solution of the initial-boundary-va¬ 
lue problem (41) - (43) 


“ t > / ii -Kk (t-t') 

G(x,t) - - % l K [f e n dtfjsin k X, -a < X < +a, t > 0 


dt 


(55) 


By Eqs. (33), (34), (38), (44), and (55), the magnetic field in the moving conductor 
is 
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B(x,t) =■ B o + i^t)B o /2(k+ac)](x-x(t)) - J Bq ? K J 

n*l 0 

x(t) - a < x < x(t) + a , t > 0 


]sink n (x-x(t)). 


(56) 


This Is a fundamental result, which shows that the transient magnetic field B(x,t) is 
the sum of a field F(x,t) Induced by the motion v(t) and a field G(x,t) induced by the 
acceleration dv(t)/dt of the conductor in the external magnetic field B q . Similar de¬ 
compositions exist for the electric field E(x,t) and current density j(x,t) in the ac¬ 
celerated conductor. 

The still unknown forms (±) of the two external wave functions *? + (t + x/c + a/c) 
in the spaces x ^ x(t) ± a are determined from the solution (56) by means of the boun¬ 
dary condition (20), which gives 

r (t 2 =* B(*»±a,t) , t > 0 . (57) 

± c 

The transformations t^ » t + (t) and its inverses t * f + {t + ) for the two (±) waves defined 
by 

t » t + , t - f{t + } , t>0 , (58) 

i c x x 

where t + * t_ * 0 for t * 0 (x(t*0) * 0) but t + i 4 t for t > 0, show that the wave func¬ 
tions are of the form (t + )-B(k**± a,t - f + {t + >). Accordingly, 

f± (t + - B(X-ta, t-f + { t + ) , *<t) ± a5xSi(a + ct) , t > 0. 

(59) 

Since B(x,t) * B(x,t) - B q , substitution of Eq. (56) into Eq. (59) gives the wave func¬ 
tions ¥ + as functionals of f + {t + x/c + a/c) : 
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f + (t + 2L ^~ ± ) - ±[aB /2(< + ac) ]v(f + (t + *-±-^}) 


+ — B J K sin k a e n ° 

~ k o “. n if 
n=l 




r* r t s‘ **,• 

/ f 1 - - *•. 


x(t) ± a S x ^ ±(a + ct) , t > 0 


. (60) 


The solutions (60) determine the propagation of the emitted electromagnetic waves 
outside of the moving conductor, x ^ x(t) ± a. Again, these waves consist each of a 
"velocity" wave [v(t)] and an "acceleration" wave [dv(t)/dt]. They satisfy all boun¬ 
dary and initial conditions, 

Vx - £(t) ± a " »(«*.«> * t > 0 , (61) 

l±/ t m o - 0 , x ^ ±a , (62) 

since t + » 0 and x = ±a [Eq. (60)] for t = 0, and, hence, f + ■ f + (0) * 0 [Eq. (58)] 
and v(f + ) ■ v(0) * 0. 

For a brief illustration of the transformations (58) consider the simple con¬ 
ductor motion v(t) * v , x(t) * v t, t > 0. In this case, t. ■ (1 + v /c)t and 
o o ± o 

V 

t * t + /(1 ^ —2), i.e., f + (t + ) is proportional to t . 
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ANALYTICAL SOLUTIONS 


For the most general representation of the electromagnetic fields In the moving 
conductor and the ambient medium, dimensionless Independent and dependent variables 
are introduced by 


£ - x/a, t - Kt/a , £(x) - x(t)/a, a " k a, v(x) * v(t)/v 

n n c 


(63) 


8(£,x) = B(x,t)/B Q , E(£,t) - E(x,t)/(<B q /a), j(£,x) - j(x,t)/(B Q / U a), 


^ <C « 'i (x,t)/B 


(64) 


1. Conductor Solutions 

According to Eq. (56), the dimensionless electromagnetic fields 8 (£,t), /(£,t) 
38(£,x)/3£, end E(£,t) in the accelerated conductor are: 

B(£,x) - 1 + 2(A" R ) 

* X j ,,/ls -of (t-T*) 

” M I K n f/ e dr' ]sina Q [E -|(x)] , 

n«l 0 


£(x) - 1 < £ < £(x) + 1 , t > 0 


, (65) 


^’ T > ‘ 2(1+ TZ) v(t) 




“ W J a K [/ e n dx'Jcosa [£ -£(x)] 

n-1 n n 0 T n 


£(x) - 1 < £ < £ (x) + 1 , x > 0 


, ( 66 ) 


E(£,x) - -Mv(x)B(C,x) + /(£,x), £(x) - 1 < £ < $(x) + 1, x > 0 , (67) 


where 
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- R , 


n = 1,2,3, ... 


( 68 ) 


and 


a ctga 
n n 


K - (R 2 + a 2 )sina /a 2 [(R 2 + a 2 ) + R] 
n n n n n 


M ■ ycrav 

o 

R ■ yoac 


, (69) 

, (70) 

. (71) 


M is known as the magnetic Reynolds number of the conductor with characteristic speed 
-v q . Eqs. (65) - (67) Indicate that M/(L + R) determines the order of the ratio B/B q 
of induced and external magnetic fields. The steady-state Induction in moving con¬ 
ductors ^’^ 2 ^ is determined only by M. 

- 1/2 

R is a new dimensionless group which involves the velocity of light c * (eu) 
so that R » M for v q « c. R has the physical meaning of a "magnetic Reynolds num¬ 
ber" of the ambient nonconducting space, i.e., R is a coupling parameter between the 
conductor (0 < o < ») and its external medium (cj ** 0), which determines the magnitudes 
of the external electromagnetic transients [Eqs. (75) - (76)]. 

The net electric current flowing through the conductor is per unit width An * 1 

C(t)+1 

I(t) ■ / j«-,T)de - Mv(t)/(1 + R) 

S(t)-1 

2 , 

» T./,\ -O(T-T) 

- 2M I sim K / e n dr’ , r > 0 . (72) 

L . n n > d t' 

n-1 0 


2. External Solutions 

For the space surrounding the moving conductor, the wave functions ^ + (x,t) of 
the electromagnetic transients are by Eqs. (60),(63) - (64) in dimensionless repre¬ 
sentation 


87 






By Eqs. (18) - (19). and (73) the dimensionless solutions for the electromagnetic field 
outside the moving conductor are: 

8 + (?,t) - 1 + * ± (t + , £(t) ± 1 * C s ±(1 + Rt), t>0 , 

-1 , 1 + Rr<±e<«, r > 0 (75) 

and 

E + (5.t) * +R<C ± (t + -5-±-i) , $(x) ± 1* 5 * ±(1+Rt) , t>0, 

« 0 ,1 + Rx<±5<® , t > 0 . (76) 
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DISCUSSION 


The magnetic Reynolds number is M * ycrav o 5 I for conductors (depending on the 
parameters a, a, and v ), whereas the coupling number R = y Q oac >> 1 is large for macro¬ 
scopic (a) conductors. In the most general case, electromagnetic induction in a mo¬ 
ving conductor is determined both by U and R, where M « R since v q << c. 

4 -1 -4 

1. Case R » 1 . Por most macroscopic conductors, it is a > 10 0 /m and a > 10 m 

—7 8 2 

(p - 4ir*10 Vsec/Am, c * 3x10 m/sec) so that R > 10 , and by Eqs. (68) and (69) 

a Q a mr , R >> nr , n - 1,2,3, ... , (77) 

K n * <-l) n+1 (nTr)" 1 R" 1 , R » *» , n -1,2,3,... , (78) 

It is seen that the "velocity" and "acceleration" fields In Eqs. (65) - (67) are of the 
same order since 


M/(l + R) - M(+K x ) - M/R = v Q /c « 1 , R » 1 . (79) 

In the magnetic field solution (65), the dominant term is the external 8 q - 1 >> M/R, 
and in the electric field solution (67) the dominant term is the motion induced field 
|Mv 8 q | >> I/I , if R >> 1. Although significant electric fields are induced in the con¬ 
ductor, the Induced magnetic field is small, 8 - M/(l + R) « B ■ 1 if R » 1, and, 
therefore, the current density / - 3B/3£ is small, too. 

2. Case R » °» . In an actual experiment, R - « can never be reached but only 
asymptotically approached. In this hypothetical situation, the conductor [Eqs. (65) - 
(67)] and external [Eqs. (75) - (76)] solutions reduce to: 


and 

8(C»t) 


Btt, t) - 1 , jtt.t) - 0, E(5,t) - -Mv(t), |(t) - 1 «S l < |(t) + 1, t > 0, 

(80) 

1 + OfMR" 1 ], E ± (5, x) - MR 0 ft 0 , C (t) ± 1 * C * ±(1+Rt) , t > 0 , 
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( 81 ) 



'<*• ■ 


since 


lima * nr , lim K * 0 , lim RK = (-l) n+ ^/mr, n = 1,2,3,... 


(82 


R-h» “ I?-*® “ R-*« 

It should be noted that for R-*• ®, only B + = 0 but E + 4 0, i.e., the external 
electric transients are (behind their wave fronts) of the same order of magnitu- 

-t- -2 

de as the electric field £ in the conductor . Since 7 xjj + * c 3E + /3t, the elec¬ 
tric transients E + always coexist with (no matter how small) magnetic transients 
8 . For this reason, the limit R * ® has no physical meaning, quite apart from 
the fact that always R < ® for o,a,c < ®. The conventional boundary conditions^’^’^ 
for electromagnetic diffusion processes in conductors imply R * 00 and E + s 0, and 
are, therefore, physically incorrect. 

As an explanation it is noted that, in conductors, magnetic field diffusion 
is a nonrelativistic process, as is electric conduction, j * o(E + vxB). The elec¬ 
tric transients E + in vacuum must be of the same order as the electric field E in the 
conductors, E + ~ E ~ MR®, since otherwise the tangential electric field would not be 
continuous across the conductor surface. On the other hand, the external magnetic 
transients 6 + are small of order MR ^ ■ M(p Q aa) ^c ^ since the magnetic field energy 
flows with the speed of light in the ambient space. The deeper physical reason for 

the external electromagnetic transients is to be seen in the conservation laws for 

9 ) 

electromagnetic energy and momentum, which follow from Maxwell's equations 

3. Diffusion Approximation . It Is known that Maxwell's equations with dis- 

placement current and the nonrelativistic Ohm’s law, j ■ a(E + v*B), combine to a 

-*■ 9 ) 

hyperbolic diffusion equation for the magnetic field B in conductors , which reads 
in the considered one-dimensional case with r-independent conductor velocity v(t) 


dh _1_ 3B v_ 3B 2 3 2 B 
3t 2 T R 3t T R 3 *" C 3 x 2 
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"W"B» 


( 83 ) 



where 


t = e /a << 1 (84) 

K. O 

is the field relaxation time, which is extremely small for conductors. Eq. (83) 
reduces to the parabolic diffusion equation in the limit t r << 1: 


3B 3B 
3t ^Tx 


2 

2 3^B 

T » c 77 


t r “ 1 


(85) 


The parabolic diffusion equation is an excellent approximation, since the relaxa¬ 
tion time of conductors is very small, r R << 1. By Eqs. (84) and (85), the field 
relaxation time -t r and the diffusion time t d are interrelated by 


-1 


2 / 2 
_c /a 


( 86 ) 


where a is the extension of the conductor. For conductors, the diffusion time is 

relatively large if a is not too small, i.e., t d >> t r . 

2 2 

Comparison of the neglected term 3 B/3t with the leading (0 < |v| << c) 
second and fourth terms of Eq. (83) reveals the relation of the parabolic diffusion 
approximation to the new coupling number R ■ u Q oac: 


|3 2 B, ,i -1 3Bi i3 2 Bi,i 2 3 2 Bi T R e o P o 

1 — 2> 7 1 t r at ' ! 7 c —7 -7"* - ? 

3t 2 R 3t 3t 2 3x 2 T D ( W aa) 2 


,-2 


(87) 


This result again confirms the validity of the parabolic diffusion eouation for 

conductors, for which R - p Q oac >> 1. More important, Eq. (87) demonstrates that 

2 2 -2 

the neglected relativistic term 3 B/3t is small of order R «<< 1, whereas the 
calculated electromagnetic fields in the conductor are of order 8 - B - MR? [Eqs. 
(65) - (67)], and the external electromagnetic transients are of order B + ~ MR ^ 
and E + - MR 0 [Eqs. (73) - (76)]. 
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Thus, consistent electromagnetic field solutions for the regions in and 
outside of the accelerated conductor have been obtained within the parabolic 
diffusion approximation, which satisfy the boundary conditions for the conti¬ 
nuity of the tangential electric and magnetic fields. If an accuracy of higher 
_2 

order R is to be achieved, then the hyperbolic diffusion equation (83) has to 
be used. For non-relativistic conductor motions, ]v] << c, however, an accuracy 
of order R * is completely sufficient. The mathematical advantages of the para¬ 
bolic diffusion equation become obvious if it is used in connection with (time- 
dependent) moving boundary conditions (magnetic flux compressors, electromagnetic 
induction generators, etc.), which are extremely difficult to treat for the hyper¬ 
bolic diffusion equation. 
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ABSTRACT 


The heating rate of electrons in laser-irradiated plasmas is derived 
from the quantum-mechanically extended Vlasov equation. The heating of 
electrons is shown to be due to (i) the electron velocity space diffusion 
by the plasma waves, which are induced by the laser field (classical tur¬ 
bulent heating), (ii) the microscopic quantum-mechanical interaction 
between each electron and the laser field in the electric field of the 
plasma waves (collisionless inverse bremsstrahlung), and (iii) collisio- 
nal inverse bremsstrahlung in the static plasma field. A quasilinear equa¬ 
tion for the changes of the electron distribution function by these three 
concomitant mechanisms is derived and the heating rates of the electrons 
are calculated. It is shown that the collisionless inverse bremsstrahlung 
is dominant for hot plasmas in most experimental cases. The so-called anomalous 
heating is identified as being due to collisionless inverse bremsstrahlung, 
and is explicitly expressed in terms of the wave number of the plasma wave 
and the plasma parameters, when the turbulence is due to the ion acoustic 


instability and is stabilized by ion trapping. 
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INTRODUCTION 


If a laser beam is injected into plasma, the energy of the laser field 
is transferred to the electrons by the following mechanisms: 

(i) In part, the laser radiation is transformed into plasma waves through 
interactions with electrons or the already existing plasma waves. Then, 
energy flows from plasma waves to electrons by Landau damping, which diffuses 
the electrons in velocity space. This mechanism is known not only as the only 
one heating mechanism caused by the plasma waves but also as the responsible 
mechanism for the anomalously large heating of collisionless plasmas by la¬ 
sers [1]. Hence, this mechanism has been called either turbulent heating 

or anomalous heating by other authors [1,2,3]. However, we call it "classi¬ 
cal turbulent heating" since this mechanism turnes out to be neither the sole 
heating mechanism of the plasma waves nor the dominant cause for the so-called 
anomalous heating. 

(ii) The direct absorption of the laser radiation by the individual electrons 
under the influence of the plasma waves (collisionless inverse bremsstrahlung) 

[4]. 

(ill) The same process as collisionless inverse bremsstrahlung, but under the 
influence of the static plasma field (collisional inverse bremsstrahlung) [5]. 

These three mechanisms are independent but can occur concomitantly as 
one readily sees from the quantum-mechanically extended Vlasov equation (Kim - 
Wilhelm equation). 

The first systematic theory of the dissipation of an electromagnetic 
wave in plasmas has been given by Dawson and Oberman [6]. However, since their 
theory is based on the hydrodynamic momentum equation (force equation) of the 
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electrons only, the heating of the electrons cannot be treated by their theory 
in principle (only the over-all ohmic dissipation of the electromagnetic wave 
energy, including scattering, can be determined). Furthermore, the application 
of their theory is limited in the weak laser-intensity regime where anomalous 
heating by laser-induced plasma waves does not exist. 

Recently, Manheimer [7] has conjectured an equation for the heating of 
the electrons in the laser-irradiated plasma by adding a phenomenological term 
for the anomalous heating to the energy equation of the electrons [7,8]. The 
latter appears to have been guessed from the Vlasov equation without taking in 
account properly the laser field. The term for the so-called anomalous heating, 
which is supposed to consider heating by collisionless inverse bremsstrahlung, 
is purely phenomenological, i.e., is not derived from first principles. 

For these reasons, herein the equation for the heating rate of the elec¬ 
trons of the plasma in the field of high-frequency electromagnetic waves is 
derived from the quantum-mechanically extended Vlasov equation. The quantum- 
mechanically extended Vlasov equation is obtained from the ordinary classical 
Vlasov equation by replacing the term for the change of the electron distri¬ 
bution due to acceleration (or deceleration) of the electrons by the laser 
radiation through a term representing the classical limit of the change in 
the electron distribution due to absorption (or emission) of laser quanta [9]. 
The equation for the heating rate of electrons here gives the detailed features 
as to how and what amount of the energy of the laser radiation is transferred to 
the electrons. Especially the dominant mechanism for the anomalous heating of 
collisionless plasmas by lasers is clearly revealed by this equation. 
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BASIC EQUATIONS 


We begin by adopting the set of Vlasov equations to describe the change 
of the electron distribution in the (r,v) space, and the time-dependent Schroe- 
dinger equation to describe the dynamics of the individual electrons in the r- 
space. The ions are regarded as a set of randomly distributed fixed scatters. 

We assume that the laser light is a circularly polarized electromagnetic 
wave propagating in the z-direction. The spatial dependence of this electro¬ 
magnetic wave is neglected since the wave length of the electromagnetic wave 
is assumed to be much longer than that of the plasma waves. Under these 
assumptions, the field equations of the plasma are (Gaussian units): 


. 21 - £ (1 - V$).— 

3t ^ m e.m. ^ 


+ v 


V (p m 4jr[e/f d v - Ze£s(r - r.)] 

j 2 


4« f£ - t (l/2m) (| V + | l) 2 - e* 


where 


ial 


e.m. c3t 


, ( 1 ) 

, ( 2 ) 

, (3) 

, (4) 


A(t) * A q (x coswt + y sinuit) . (5) 

4>(r,t) is the potential of the plasma field (the plasma waves plus the static 
plasma field). The subscript e.m. is an abbreviation for "electromagnetic", 
and the other symbols have the usual meaning. 






The potential of the plasma field can be expanded in a double Fourier 


series, 

4>(r,t) - l £*<£,fl)e 1( ^' r - nt) . (6) 

l n 

(For the case that the static plasma field is negligible, ij>(r,t) can be ex¬ 
panded in a single Fourier series since in this case tc and ft are interrelated 
by the dispersion relation). Since <t>(r,t) is real, we have <p(—tc, —ft) * <p(tc»ft)*. 

After substituting Eqs. (5) and (6) into Eq. (3), we solve for the tran¬ 
sition probability amplitude from the resulting equation by means of first- 
order perturbation theory [4,5] where ip is the pertubing potential. The tran¬ 
sition probability amplitude from a state 1 with momentum p^ to a state 2 with 
momentum p 2 of the electron is [4,5] 


a(l 2) 


« . . sin[ (— - ft - nw)T] 

■r « V>'“> ■ < 7 > 

£ ft n—<» 


n#0 


where 


X * eE Q k i /mu , e - (p^ - p*)/2m . (8) 

6 is the azimuthal angle of i.e., k * (k^,6,k z ), and c' is a real constant, 
which drops out since exp(-ic') does not affect the transition probability. 

The subscript i refers to the direction perpendicular to the direction of pro¬ 
pagation of the electromagnetic wave, and E q is the electric field amplitude 
of the electromagnetic wave. 

From Eq. (7) we see that the transition probability per unit time for the 
transition from state 1 to state 2 (with the absorption or the emission of pho¬ 
tons) is 
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(9) 


T(1 - 2) =■ T(v lt v 2 ) = |a(l - 2)| 2 /2T 


- ^ X I ^ I 2 I jJ(A/u)«(e -Itt - nftu) 

t n n= -°° 

* n yo 

In the derivation of Eq. (9), we assumed that there are no phase relations be¬ 
tween different components of the plasma field. 

The change of the electron distribution due to the absorption or emission 
of photons of the electromagnetic wave [4,5] is 


[3f(v )/3t] • 

i e .m. 


Ira 

i 


2)f(v x ) - 7(2 


l)f(v 2 )] 


Fd 

1 


1 ,v 2 )[f(v 1 ) - *(v 2 )] 


where 


( 10 ) 


v, * p./m ■ 2'ir'tim ^(n /L ,n /L ,n /L ) , n ,n ,n * integers, 

11 x x y y z z x y z ° 

is quantized (L x ,L y ,L z are the x-,y-,z-dimensions of the plasma, respectively). 

In Eq. (1), (e/m)^ -3f/3v is the rate of change in the number of electrons 

e. m. 

per unit volume in the (r,v)space due to acceleration (or deceleration) of elec¬ 
trons by the laser field. From the quantum-mechanical viewpoint, the acceler¬ 
ation of an electron by the spatially uniform laser field is possible only by 
absorption (or emission) of the laser photons in the presence of another field 

which is the plasma field in our situation. Hence, (e/m)l *3f/3v is equal 

e .m 

to the classical limit (h -*• 0) of the rate of change in the mumber of electrons 
per unit volume in the (r,v) space due to the absorption (or emission) of the 
laser photons by the electrons under the influence of the plasma field ("inverse 
bremsstrahlung"). Accordingly, by using Eqs. (1), (3) - (6), and the same meth¬ 
ods as in Ref. 5 and 6, we obtain the quantum-mechanical Vlasov equation, which 
is valid up to second order in $: 
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(11) 



fl+ntO v 3 f (y) 
U ' o 


], 


3v_ 


t 

( 12 ) 


with * v*1c/k. The subscript i.b. is an abbreviation for "inverse bremsstrah- 
lung'.'. In the derivation of Eq. (12), v + * -v >* and <f>(“k»“0) * (j>(tc,n)* have 
been used to cancel out the terms, which change their signs as (1c,fi,n) changes 
to (-1c,-fi,-n), in the summation over tc.fi, and n. 

Equation (11) is the desired Vlasov equation of the electrons in the 
electromagnetic field and the plasma field, in which the term (3f/3t)^ con¬ 
siders all interactions of the electrons with electromagnetic field. It should 
be noted that the presented derivation of Eq. (11) is based on first principles 
(Vlasov and Schroedinger equations in electromagnetic and plasma fields). 

Since Eq. (11) is valid up to the second order in ij>, we can apply to this 
equation an analogy to the derivation of the quasllinear equation for the electron 
distribution in the absence of the electromagnetic wave [10,11]. The result, 
i.e., the quasllinear equation for changes in the electron distribution in the 
presence of the electromagnetic wave, is by Eq. (11) 


3f 

at 


( —, + (ii) 

v 3t'c.t. 


at'i.b. 


(13) 


where 


(—) 

v 3t'c.t. 


- 41 

* 


_3 .[ lc(it-3f/3v) j 
3v - fi - iy 


(14) 


100 










is the change in the electron distribution associated with the diffusion in 
the electron-velocity space due to the acceleration of the electrons by the 
electric field of the plasma waves (or the change in the electron distribu¬ 
tion by classical turbulent heating). The subscript c.t. is an abbreviation 
for "classical turbulent heating", and y > 0 is infinitesimal. In Eqs. (13) 
and (14), the spatial average of the actual distribution function is desig¬ 
nated as f, and the same simplification in notation is adopted in the fol¬ 
lowing, where there is no chance of confusion. 
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HEATING RATE OF ELECTRONS IN COLLISIONLESS PLASMA 

The heating rate of the electrons per unit volume under the influence 
of the external electromagnetic and internal plasma fields is by microscopic 
definition 


3<e> 

at 


,3 1„ /* »2 r , 3fv 

d - u) [(^) 1#b> 


St'c.t. J 


( 1 ^) 

v at 'c.t. 


(15) 


where N is the unperturbed electron density. In Eq. (15), the first term is 
the heating rate of the electrons per unit volume by inverse bremsstrahlung, 
and the second term is that by classical turbulent heating, i.e., the rate of 
energy transfer to the electrons from the plasma waves by velocity-space dif¬ 
fusion. 

We model the spatial average of the electron distribution as 


f(v) - (2,„V 3 'W-k^l{l ♦ ( t£# L 2 - 3)} (16) 

2v { 3Nmv v (r ' 


where the hydrodynamic velocity u and heat flux 0 are induced by the laser field. 
In the appendix we calculate the second term of Eq. (15) as 


(2f-i c , - <f> 1/2 *4 1 i lt( ¥ )j2 - t-s> 

J * mv f n 2v k 

e k e 

f.. ■* -*• . (£•$)■* „,n^£*u .2 . 1,(1 - ic'Ux4A 

T " *' u 1 " ( “^— 5 3 ( -w—> 1 * (17) 

' Nmv k Q “ “ 

e 


The calculation of the first term by means of Eqs. (12), (15), and (16) 
is elementary, and gives 
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• -4i i 2 /i ££•«<* - : > 2 >“ 

ra t 0 J n*-® K t t 


k k 

(7) 1/2 l J J 2 (i ) exp[--^f|^ ] (Q-Hw - bZ) 

2 3 L k k L nu r » 2, 2 

mv^ £ 0 n*-<= 2v k 


nfO 

'id 2 

k u n 'oj 

n=-<x 

ni<0 

, 3 


L * m - i-s + u . 2 ( 2jssi»,2 + . (18) 

' Nmv 


Equations (17) and (18) can be used irrespective of the origin of the 
plasma field as long as the strength of this field is so small that the resul¬ 
ting perturbation of the electron distribution can be adequately treated by 
means of first-order perturbation theory. Bearing this in mind, we reduce 
the expression of the heating rate to a simpler form for the cases of small 
and large A/w. 

(i) A/oi << 1 ; In this case, the laser-induced current and heat flux 
can be neglected. Then, by means of the same approach used to derive Eqs. (1) 

- (15) of Ref. 6, we ontain from Eqs. (1) and (2) 


. ,■* .x (2it) 3 v i£*r r „/?-x „,+x* -i6 iwt . i<5 -iuit, 

<|>(r,t) « -*—le [S(k) - C(k)"e e + C(k)e e ] 


where 


with 


S(t) 


Ze 


2* 2 (k 2 + kj) j 
-ilt’ri 


—±tc*r. 

h 3 


C(tr) m -5s_ A. (Vo j-v k ~ 2 r_i_i_I* 

() 4w 2 u Y )k D(k,0) D(k,u) 


D(k,«) - 1 + -yW(s) , s - -f-r 
k e 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 
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Here, D(k,w) is the longitudinal dielectric constant, V is the volume of the 


plasma system, = oi^/v^ c ^ e De ^y e wave number, the plasma frequency, 
and W(s) is the plasma function given in the appendix. 

Designating by $^(r) the Coulomb field potential at r due to an ion at 
rj, and by X Q the Debye length, the first term of Eq. (19) becomes 


<Kr) - I<Mr) =• l exp(-[r-r | A ) , (23) 

j JI r ~ r jI J 

Hence, the first term of Eq. (19) is the static Coulomb field potential, re¬ 
duced by the usual Debye shielding, acting on an electron due to all ions in 
the plasma. This static plasma field gives rise to collisional inverse brems- 
strahlung but does not contribute to classical turbulent heating. The second 
and third terms are the potential of the laser-induced plasma waves, which can 
contribute both to collisionless inverse brensstrahlung and classical turbulent 
heating. 

From Eqs. (18) - (21) and 

J 2 (~) - —^-(T-) 2 ^ for X/oj « 1 , (24) 

n a) f | v 2 2 to 
(n!) 

it is seen that the collisionless inverse bremsstrahlung is of the fourth 
order in X/w, and can be neglected compared with the collisional inverse brems¬ 
strahlung, which is of the second order. Hence, up to the lowest order in x/u>» 
we obtain from Eqs. (18) - (21), and (24), 


:e>v . f „ x 7/2 Ne 2 fd 3 k 1 |S(£)| 2 X,2 2 , . 2 ,-2.2 

rr\.b. (2,r> ~l~ kY ^ u exp(_u /2v e k > 

mv J 
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The expression for tpj(r) given in Eq. (23) shows that <pj(r) decreases 
very fast to zero as | "r-r.. | /X Q ». This decrease is mathematically brought 
about by representing discrete phenomena through smoothly continuous functions; 
<t>j(r) is in reality zero for |r-Tj| >> Xq. This corresponds to S(k) * 0 for 
k << kp. Considering this behavior of S(tc) and Eq. (20), we can explicitly 
evaluate Eq. (25) as 

„ . ,,, ,2 6.2„2 ^ k D dk k 3 exp(-w 2 /2v 2 k 2 ) . . Z 2 e 6 E 2 N 2 

, 3<e>> - A/« >1/2 Z e E„N / ___ e__ _ 4/o >1/2 o T! ,„ > 

( at \.b. 3 <2 5 3 3 2 y ..2 , ,2.2 ~ 3 (2ir) 3 3 2 U(z) ’ 

m v e “ 0 < k + V “ V e“ 


U(z) ■ / dx[x 3 exp(-l/2x 2 )/(x 2 +z 2 ) 2 ] , z * w /w > 

0 P 

with ojp the plasma frequency. The function U(z) is represented in Fig. 1. 

In the derivation of Eq. (26), we assumed the ion distribution to be 
random so that 


^1 l exp(-ik-Tj)j 2 » N 


since the large-amplitude ion wave cannot be induced by weak laser fields for 
which X/u << 1 [12]. 

From Eqs. (17) and (19), we obtain for the heating rate per unit volume 
by classical turbulent heating, up to the lowest order in X/u, 


,3<e>> - >7/2 e 2 N rd\[c(t )| 2 2 , 2,„..2, 2, 

( nr ) c.t. (2,t) —3 /— V 10 ex P (_w /2v e k > 


• (29) 


Eq. (29) can be expressed more explicitly for the cases of hot and cold plas¬ 
mas. In this connection, it is noted that electromagnetic waves propagate 
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Fig. Is U(z) versus z by Eq. (27) 







in plasmas only when their frequency is greater than the plasma frequency. 
Accordingly, 

v g < wv e /u) p = u/kp . (30) 

For hot plasmas for which >> w/k, Eq. (22) becomes, to the lowest or¬ 
der in w/kv^, such that [1/D(k,0) - 1/D(k,ai)] is not zero, 

D(k,oj) = 1 + k^/k 2 - i(ir/2) 1/2 (k^/k 2 )^ . (31) 

e 

Also, from v^ >> u/k and Eq. (30), we have k >> k^. This means that the 
plasma behaves, as far as the plasma waves are concerned, like a system of free 
individual particles [13]. 

Substituting Eqs. (21) and (31) into Eq. (29) and using k >> kp yields 


Z 2 e 6 E 2 N 2 


/ a< e > N , JL,- .1/2 u e l o N f dx x 

v 3t c.t. 3 V ’ 3 3 2 2.4 

m v a) ~1 (1 + x ) 
e p 


, 0 ,3/22 622 
(2tt) Z e E N 
o 

3 3 2 
144m v w 
e p 


. (32) 


By comparison of Eq. (32) with Eq. (26), it is recongized that the hea¬ 
ting rate of the electrons by inverse bremsstrahlung is slightly greater than 
that by classical turbulent heating for &>■<«» . 

For cold plasmas for which v fi < w/k^ << us/k, Eq. (22) becomes, to the 
lowest order in v e k/w such that [1/D(k,0) - l/D(k,co)] in Eq. (21) is not zero, 

D(k,u)) ■ 1 - w 2 /u) 2 . (3: 


Accordingly, Eq. (29) reduces to 


3t c.t. 3 


T < 2*) 1/2 


Z 2 e 6 E 2 N 2 w 4 


mW (w 2 - u 2 ) 2 
e p' 








where 


V(z) - / dx[exp(-l/2x )/x(x 2 +z ) ] . (35) 

0 

The function V(z) is represented in Fig. 2. 

Equation (34) appears to indicate that classical turbulent heating by 
the plasma wave becomes resonant as u + However, Eq. (34) should be ques¬ 

tioned since Eq. (38) is derived by using Eqs. (19) and (21), which are derived 
from the collisionless Boltzmann equation (or Vlasov equation), even though the 
considered plasmas are cold, and consequently, collisional. 

Comparison of the expression for the heating rate by inverse bremsstrah- 
lung, i.e., Eq. (26) with that found by Seely and Harris [Eq. (15) of Ref. 5] 
shows that our result is the same as theirs except for a factor U(z)/lnA, where 
InA is the Coulomb logarithm, 

InA - ln(12irNAp) . (36) 

Since the Coulomb logarithm for laboratory plasmas is about 10 to 20 and U(z) 
is less than 0.04, the value of the heating rate given by Eq. (26) is less than 
about four-thousandth to two-thousandth of that given by Eq. (15) of Ref. 5. 

This difference is due to the fact that Seely and Harris do not consider Debye 
shielding, whereas our result contains Debye shielding brought about by the 

collective dynamics of the electrons [Eq. (19)]. 

2 

(ii) A/w ■ eE^/mu >> 1: For this case, we can use no longer Eqs. (19) 
- (21) for the plasma field potential $ even if there is no other perturbation 
except the electromagnetic field. However, we can readily surmise that the 
collisio nless inverse bremsstrahlung is stronger than the collisional one for 
y s A/w » 1. This is because the second and third terms of Eq. (19), which 
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give rise to collisionless inverse bremsstrahlung, become important compared 
to the first term, which produces collisional inverse bremsstrahlung, as y 
increases (but still y < 1). For this reason, the result obtained by Seely 
and Harris for the effective collision frequency in the case y >> 1 (Eq. (22) 
of Ref. 5] is questionable since they take into account only collisional in¬ 
verse bremsstrahlung. 

Except for increasing y (usually by increasing the intensity of the elec¬ 
tromagnetic field), the collisionless inverse bremsstrahlung becomes more impor¬ 
tant compared with the collisional one with increasing electron temperature. 

This is because the mean free path of the electron becomes larger with increasing 
electron energy. Accordingly, the collective effects which are enhanced by a 
large number of instabilities (not caused by the local electromagnetic field) 
become more dominant than the effect caused by the static ions. Hence, even if 
y < 1, the collisionless inverse bremsstrahlung may dominate the collisional one 
for very hot plasmas if there are causes of turbulence other than the local elec¬ 
tromagnetic field. The latter case will be treated in the next section. 

For y * X/w >> 1, we can write approximately [5] 


u L J l&- Mr 4 * ■ T'« <* t - “-M* + «<» t - MM 1 • <37 > 

n^O 

Substituting Eq. (37) into the first relation of Eq. (18), we find for hot 
plasmas, v g >> Q/k, X/k, 


<Tf>i.b. * <f> 1/2 *4 I {to - t.S)[a - t-Z * -&&L] + 444} . 


Nmv k Q mm* 











Comparison of Eq. (38) with Eq. (17) indicates that the heating of the 
electrons by collisionless inverse bremsstrahlung is greater than that by 
classical turbulent heating where the following conditions hold: u >> 1 and 
v g >> fi/k, X/k (hot plasmas). 
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ANOMALOUS HEATING RATE 


Heating by coillsionless inverse bremsstrahlung, which is the main heating 
mechanism in laser fusion, is usually called "anomalous" heating. This phenomenon 
requires a correct quantum-mechanical treatment of collisionless inverse brems- 
strahlung, as distinguished from the deficient classical approaches in Refs. [6,7, 
8,12,13 ]. From the experimental point of view, the heating rate is "anomalous" in 
collisionless plasmas since it is unexpectively large compared with the classical 
theories, which consider collision-dominated phenomena. In the following, the 
significance of coillsionless inverse bremsstrahlung is demonstrated by comparison 
with the classical theory. 

Equation (38) is valid irrespective of the causes of the electrostatic field. 
In laser fusion, it is believed that the electrostatic field is developed by the 
ion acoustic instability driven by the return current [7], For this case, it is 
known from computer simulations that the R-spectrum of the electrostatic field is 
flared out in a cone of 45° to the direction of the return current [7,15-17]. Hence, 

we can assume that is of the order of k. Then, the condition p > 1 corresponds 

12 2 

to a laser intensity of I > 10 watt/cm for a neodymium laser. With presently 
available laser powers, this condition is generally satisfied in most experiments 
so that Eq. (38) is applicable to them. 

Both computer simulations and analytical theories indicate that the ion 
acoustic instability is "stabilised" by ion trapping [18-22]. Hence, its intensity 
level can be estimated on the basis of simple trapping arguments to be [7] 


KT 

TT- t 


3T 


<k 2 ; kj) i/z 


v 


(39) 


where the wave spectrum of the turbulent electric field is modeled for a single 
wave having wave number k, frequency 0, and a propagation angle of 135° with the 






direction of the laser beam. By incorporating the estimate (39) for the 
saturated turbulence level into Eq. (38), we obtain the anomalous heating 
rate for hot plasmas with v >> A/k: 



where Np is the number of electrons in a Debye sphere [23]. By definition, in 
"collisionless plasmas" the collective behavior far dominates the collisional 
effects [23], and N Q >>> 1. Thus, by Eq. (41) we demonstrate the significance 
of our collisionless quantum-mechanical theory compared with the classical col¬ 
lisional theories and explain the observed anomalous heating rate in experiments, 
R i 6 N d »> 1. 
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CONCLUSION 


Under che condicion that the field strength of the electromagnetic wave is 
below the threshold above which strong turbulence is developed (so that the first- 
order perturbation theory would no longer be appropriate) we have obtained the fol¬ 
lowing results: 

(1) From the set of Vlasov and Schroedinger equations, the heating rate of the 

electrons in the electromagnetic field is systematically derived in terms of the elec¬ 
trostatic field potential, which is produced by the static ions and the collective 
dynamics of electrons in the self-consistent manner. The Vlasov equation is genera¬ 
lized by formulating in it —1 •3f/3v as a quantum-mechanical interaction integral. 

m e.m. 

(2) Especially for ji < 1, the electrostatic field potential $ in the heating 
rates is expressed in terms of the field strength E q , frequency u , and wave vector k 
of the electromagnetic wave for both hot and cold plasmas. 

(3) For very hot plasmas with v^» fl/k,*/k, it is shown that the heating of the 
electrons by inverse bremsstrahlung is considerably stronger than that by Landau 
damping of electrostatic waves. Furthermore, it is shown that for U < 1 the above 
statement is true regardless of the magnitude of the electron temperature. 

(A) The anomalous heating observed in the experiments is identified as the 
heating by inverse collisionless bremsstrahlung, and explicitly formulated when the 
turbulence is developed by the ion acoustic instability and "stabilized" by ion 
trapping. 








APPENDIX 


Combination of Eqs. (14) - (16) yields 


2/~\“l/2 


&-> c t - irl zh*^)| 2 (2-) 


</ dx ex P<~ x V 2 ) [x 2 + -i-—. &&) 3 <x 5 _ 6x 3 + 3x) ] . (Al) 

x “ z " 1Y 3Nmv J Q 


Integration over x in Eq. (Al) is accomplished by using the following 
formulae: 


(2n )* 1/2 I* ” 6 * X _ e ^- - X - -/ - 2 1 - W(z) , (A2) 

—ao X Z Y 

(2tt) — 1/2 /* " d « x 2 exp , C-x 2 /2) . zW(z) , ( A3 ) 

«ao X Z Y 

(2ir )“ 1/2 / dx x 3 e 3 g<~ lL- / 2 ) » i + z 2 W(z) , (A4) 

iaoX-Z-iy 

(2tt ) -1/2 / *x^v&iz]LJ2L , 3 + z 2 [l + z 2 W(z)] , (A5) 

' x - z ~ iy 


where W(z) is the plasma function [24] given by 
2 


i£) 1/2 z exp(-V) + [1 - z 


,4 , .,>n+l n+2 

2 + £_ .... + (-1? « _ 

3 (2n+l)M 


W(z) - 


. .x. 1/2 . z 2 . ,1.3. . (2n-l)!! . 

i(j) z exp(- —) - t— + — + * * * + 2n ' + ] ’ 


••• ],for 1z| < 1 

(A6) 

for |z| > 1. 


Substitution of Eqs. (A2) - (A6) into Eq. (Al) and obliteration of all odd 
terms for both k and ft in the summation over k and ft gives Eq. (17). 
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